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1 Introduction 

Among nonlinear partial differential equations (PDEs) there is an exceptional 
class of so-called integrable equations. Integrable equations can be studied 
with the same completeness as linear PDEs, at least in principle. They 
possess a rich set of exact solutions and many hidden properties such as 
infinite hierarchies of symmetries, conservation laws, etc. 

There are two known classes of integrable equations. One class is linearis- 
able equations, i.e. equations related to linear ones by differential substitu- 
tions. For example, the famous Burgers equation 

ut = Uxx + 2mm^. 

can be linearised by a differential substitution u = 4'x/4' (the Cole-Hopf 
transformation). In the new variable (p it takes the form of a linear heat 
equation: 

Another class is equations solvable by the inverse scattering transform method 
(such as the Korteweg de-Vries equation Ut = Uxxx + QuUx). There is mas- 
sive literature on integrable equations, their solutions and properties (see for 
example monographs [U El El HI [5] and references therein) . 
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In the symmetry approach the existence of higher symmetries of a PDE, 
or more precisely, of an infinite hierarchy of higher symmetries is regarded as 
the definition of its integrabihty. Existence of a finite number of symmetries 
of a partial differential equation may not secure its integrabihty [HI [71 El [9] . 

The symmetry approach based on a concept of formal recursion opera- 
tor has been formulated and developed in works of Shabat and co-authors 
(see for example review papers [TUl [TH [T21 [IS])- It has been shown that the 
existence of an infinite hierarchy of symmetries or local conservation laws 
or a possibility to linearise a certain equation imply the existence of a for- 
mal recursion operator. Formal recursion operator is carrying information 
about integrabihty and is not sensitive to lacunas in the infinite hierarchy of 
symmetries or conservation laws. The conditions of its existence give inte- 
grabihty conditions for the equation which can be formulated in a form of an 
infinite sequence of canonical densities and encode many hidden properties 
of the equation. If a density is non-trivial (i.e. is not a total derivative) it 
provides a density of a local conservation law of the equation. For linearis- 
able equations all densities except a finite number are trivial. The sequence 
of canonical densities is invariant with respect to invertible (and almost in- 
vertible [Hj) transformations. 

It was the first approach enabling to give a complete description of inte- 
grable evolutionary equations of the form 

Ui; f i^Uxxj Uy t) , 

ut = d^u + f{d2~^u, . . . , Ux, u) , n = 3,4, 5, 

where / is a smooth function of its arguments [151 [13 [IZ]- Integrable 
differential-difference equation of the form 

Un,t = f{Un-l,Un,Un+l) , H^Z 

have been classified in [181 [19]. It can also be applied to systems of equations. 
In particularly a complete classification of systems of two equations of the 
form 

ut = A{u)uxx + F{u,Ux), det(A(u))^0, u={u,vf (1) 

possessing infinite hierarchies of local conservation laws has been given in 
[201 [2T1 [Hj . Existence of higher conservation laws immediately implies that 
the trace of ^(u) vanishes. System ([1]) can still possess an infinite hierarchy 
of symmetries even if the trace of A{u) is not equal to zero. In the latter 
case the system is linearisable by a Cole-Hopf type differential substitution. 
In the case of A{u) being a unit matrix integrable systems of the type ([T]) 
have been listed in [22] . A complete classification of polynomial homogeneous 
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integrable system of the form ([T]) in the case when matrix A{u) is constant 
and has two distinct eigenvalues is given in [23]. 

In this approach it is not assumed that equations are polynomial or ratio- 
nal functions of independent variables and their derivatives. A disadvantage 
is that for every fixed order of the equation the integrability conditions have 
to be derived from scratch and thus it is difficult to draw a global picture, i.e., 
in all orders. Also, this approach is heavily based on the concept of locality 
which makes it difficult to apply to integro-differential, non-evolutionary and 
multi-dimensional systems. 

In this article we would like to give a brief account of recent development 
of the symmetry approach. The progress has been achieved mainly due to 
a symbolic representation of the ring of differential polynomials which en- 
able us to use powerful results from algebraic geometry and number theory. 
Symbolic representation (an abbreviated form of the Fourier transformation) 
has been originally applied to the theory of integrable equations by Gel'fand 
and Dikii [21]. Symmetry approach in symbolic representation has been 
formulated and developed to the problem of the global classification of in- 
tegrable evolutionary equations in [2Sl 1251 [2Z]. In symbolic representation 
the existence of infinite hierarchy of symmetries is linked with factorisation 
properties of an infinite sequence of multi-variable polynomials. Symbolic 
representation is a suitable tool to study integrability of noncommutative 
[2E], non-evolutionary [2niEniEI], non-local (integro-differential) [32], multi- 
component [231 [331 [31] and multi-dimensional equations [33] . It is convenient 
for testing integrability of a given system, provides useful information on the 
structure of the symmetry hierarchy and suitable for global classification of 
integrable equations. In this framework it is natural to define approximate 
symmetries and approximate integrability. Study and classification of ap- 
proximately integrable equations is a new and unexplored area of research 
with a considerable potential for applications. Symmetry approach in sym- 
bolic representation has certain drawbacks due to a restriction to the ring 
of differential polynomials, which can be amended in some cases by suitable 
extensions of the ring. 

In literature one can find other attempts to describe integrable systems 
based on properties of solutions, such as the Painleve property in the analyt- 
ical theory (see for example [21 [31 H]), existence of three soliton solutions for 
bi-linear (Hirota) representation (see the article by Hietarinta in this book), 
elasticity of soliton collisions in numerical experiments, etc. Existence of 
one nontrivial symmetry of an prior fixed order can be used for isolation of 
integrable equations [351132]. Every method has certain advantages and dis- 
advantages. The reader can judge the power of each approach by the results 
obtained, their completeness and generality. 



3 



Our article is organised as following. In the next Section [2] we give basic 
definitions and notations. We define symmetries, approximate symmetries 
and formal recursion operator. Then we introduce symbolic representation 
of the ring of differential polynomials. A generalisation to several dependent 
variables is given in the last part of this section. In Section [3l using ap- 
proximate symmetries in symbolic representation, we study the structure of 
the Lie algebra of symmetries. In symbolic representation the existence of 
approximate symmetries can be reformulated in terms of factorisation prop- 
erties of polynomials. Existence of one nontrivial symmetry enables us to 
constructively extend any approximate symmetry of degree 3 to any degree 
(Theorem H]). In symbolic representation conditions for the existence of for- 
mal recursion operator lead to a simple test for integrability. Here a state 
of art result is a global classification of integrable homogeneous evolutionary 
equations (Section 13.21) . The result of the classification can be accounted as 
following: Integrable equations are symmetries (members of infinite hierar- 
cliies) of nonlinear PDEs of orders 2, 3 or 5. Thus it is sufficient to classify 
integrable equations of order 2, 3 or 5. There is only a finite number of 
such equations (namely 10) and the corresponding hierarchies of symmetries 
(Theorem [9]). In Section H] we apply our method to non-local (Benjamin- 
Ono and Camassa-Holm type) equations. It requires a non-local extension 
of the ring of differential polynomials and symbolic representation proved 
to be a suitable language to tackle the problem. In Section [3] we present 
results on classification of Boussinesq type equations. In the case of even 
order equations we use conditions following from the existence of formal re- 
cursion operator. Together with classification results for orders 4 and 6 we 
present three new integrable equations of order 10. A global classification 
of all integrable odd order Boussinesq type equation is given in section 15. 2[ 
A generalisation of the symbolic approach to 2 -|- 1-dimensional equations 
enables us (Section [6]) to study the structure of symmetries of Lax integrable 
equations and to prove the conjecture on the structure of non-local terms. 
Finally we discuss a progress in the problem of classification of systems of 
integrable equations. In particularly we have found two integrable systems 
of order 5 which we believe are new. 
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2 Symmetries and formal recursion operators 
in symbolic representation 

2.1 Differential polynomials 

We shall adopt the following notations: Un denotes n-th derivative d''^u of the 
dependent variable u. In particular, uq denote the function u itself (often we 
shall omit the zero index of Uq and simply write u). 
A li-monomial is a finite product of the form 

where all exponents cto; are non-negative integers (a^ e Z>o); and the 

total degree is |q;| = ao + ai + • • • > 0. 

A finite sequence a — (ao, . . . , ccfe) can be seen as an element of a semi- 
group Z^Q* of infinite sequences of non-negative integers, such that only a 
finite number of entries in a sequence are non-zero and there is at least 
one nonzero entry. There is an obvious bijection between the set of all u- 
monomials and Z^q*. We can simplify the notations as follows: 

U = Uq^Ui ■ ■ ■ UjJ'. 

A differential polynomial f in variables uq, Ui, . . . with coefficients in C is 
a finite linear combination of w-monomials, i.e., 

/ = ^ aau"' , e C, 

aeA 

where the sum is over a finite set A = {a\a E Z^q*}. The set of all such 
differential polynomials is denoted TZ. It is a ring with usual addition and 
multiplication of polynomials. Moreover, it is a differential ring and C ^ TZ. 
The linear operator 

fc>0 ^ '^^ 

is a derivation of TZ corresponding to the total x-derivative. Operator 
acting on an element f & TZ results in a finite sum depending on the choice 
of /. Therefore, we do not indicate the upper limit for the summation in the 
definition of and all other operators defined in this article. 

It is easy to verify that monomials u" are eigenvectors of the following 
commuting linear operators 
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with Du{u'^) = |q;|m", and = daU"", where da = ^fc>i(^ctfc) • 

Thus the ring TZ is graded and is a direct sum of eigenspaces 

neN n,peN 

Tr = {fen I D^u) = nf} , 7^; = {/ G 7^" | x„(/) = pf}. 

\i f e n^, g e 7^™, then f ■ g e 7^;+™. Simply speaking, TZ^ is a hnear 
subspace of homogeneous differential polynomials such that each monomial 
has: (i). the number of u and its derivatives being n; (ii). the total number 
of derivatives being p. For example, 

u\u'j + 2U2U3U4 — Mq-uq G TZl C TZ^. 

In some applications it is convenient to introduce weighted homogeneous 
polynomials. Let us assume that dependent variable u has a weight A which 
is a fixed rational number. We define a linear differential operator 

Differential monomials are eigenvectors of W\ and the spectrum of Wx is a 
set S\ = {nX + m — 1 \ n, m E N} . We can decompose 7^ in a direct sum of 
eigenspaces 

Elements of are called A-homogcncous differential polynomials of weight 
II. For example, M3 + 6uui is a 2-homogeneous differential polynomial of 
weight 5. We have WnWu C W^+iy. Moreover, if A > then subspaces 
are finite dimensional. 

It is useful to define the little "o/i" order symbol. 

Definition 1 Let f e 7^. We say that f = o(7^") if f e 0fc>„7^'=. 

For example, / = o{TZ^) means that the differential polynomial / does 
not have linear, quadratic and cubic terms in u and its derivatives. 
For any two elements f,geTZwe define a Lie bracket 

[/, 9] = M9)-9M), (4) 

where the Frechet derivative for any element h E TZ is defined as a linear 
differential operator of the form 

k>0 " 
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We say that an element h G 71 has order n if the corresponding differential 
operator h^, is of order n. 

Thus 71, treated as a linear space over C, together with the Lie bracket (jl]) 
is an infinite dimensional Lie algebra of differential polynomials over C. The 
bilinearity and skew-symmetry of the bracket (jll) are obvious. The Jacobi 
identity can be easily verified. 

The grading of 71 induces the grading of the Lie algebra of differential 
polynomials since we have 

[7^^,7^™]c7^;+™-^ (6) 

2.2 Symmetries, approximate symmetries and formal 
recursion operator 

In this section, for the sake of simplicity, we give definitions suitable for 
evolutionary equations 

Ut = F, Fen. (7) 

These definitions will be later extended to non-evolutionary equations and 
the multi-component systems of evolutionary equations. 

Evolutionary equation defines a derivation Dp '-71^ 7Z: 

Dpia) = a,{F). 

In this notation the derivative = Du^ and it is in agreement with ([3]). 
Sometimes, for simplification of notations, we will denote Dp as Dt- 

Definition 2 A differential polynomial G E 71 is said to be a symmetry (a 
generator of an infinitesimal symmetry) for an evolutionary partial differen- 
tial equation ^ if the Lie bracket of F and G vanishes, i.e., [F, G] = . 

If G is a symmetry, then evolutionary equation Ur = G is compatible 
with ([7]). There are many other equivalent definitions of symmetry (see for 
example [TUl [12] )• Elements of 71 do not depend on x,t explicitly. Thus 
our definition does not include space and time dependent symmetries such 
as dilatation and Galilean symmetries. In this article, when we are talking 
about symmetries we do mean space and time independent symmetries. 

Symmetries form a subalgebra Gti{F) = {G E 7l\[F, G] = 0} which is 
the centraliser of F (it immediately follows from the Jacobi identity). Since 
F E 71 does not depend on x, t explicitly, then any equation ([7]) possesses 
trivial symmetries ui,F E Gti{F) corresponding to translations in space and 
time. 
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Definition 3 Equation ^ is said to be integrable if its algebra of symme- 
tries Cti{F) is infinite dimensional. 

For nonlinear equations of tlie form 

Ut = Un + /(Wn-l, . . . , m) 72 > 2, (8) 

it is easy to show that the algebra of symmetries C'ji{ut) is commutative. 
Moreover, the symmetries G must have a linear term (see section [3TT|) . 

Having in mind the gradation ([6]) of the Lie algebra of differential poly- 
nomials we represent the right hand side of equation ([7]) and infinitesimal 
generators of symmetries in the form 

^ = /i + /2 + ■ ■ ■ , G = ^1 + ^2 + ■ ■ • , fk,9k^n^ 
and study them in the sequence of terms, i.e., linear, quadratic, cubic, etc. 

Definition 4 A differential polynomial G E TZ is said to be an approxi- 
mate symmetry of degree p for evolutionary partial differential equation ([^ 
G] = o(7^^). 

Equation ([8]) possesses infinitely many approximate symmetries of degree 
1. An equation may possess approximate symmetries of degree 2, but fail 
to possess approximate symmetries of degree 3. An integrable equation pos- 
sesses infinitely many approximate symmetries of any degree. In the next 
section, using symbolic representation we formulate the necessary and suf- 
ficient conditions for the existence of approximate symmetries of arbitrary 
degree. 

For example, equation 

Ut = + 5uui , (9) 
has approximate symmetry of degree 2 with a generator 

G = Ut + Tuus + 14:UiU2 ■ 

Indeed, 

[G, Ms + 5uui] = 210ulu2 + 105uul + lOduuiU^ = o{lZ'^). 

Moreover, equation has infinitely many approximate symmetries of de- 
gree 2 (this fact will become obvious in the next section), but fails to have 
approximate symmetries of degree 3 and thus it is not integrable. 
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It follows from the Jacobi identity and that approximate symmetries 
of degree n form a subalgebra of TZ which we denote C^{F). Obviously 

n = CUF) D CliF) D CiiF) D ■ ■ ■ D C^{F) = Cn{F) . 

Formal pseudo-differential series, which for simplicity we shall call formal 
series, are defined as 

A = amD"^ + a.rr.^iD'^'^ + ■ ■ ■ + ao + a-iD-^ + ■ ■ ■ at ell. (10) 
The product of two formal series is defined by 

aDl o = aibD'^^^ + ClbiD^^""-^ + Clh2Dl+'^-'' + ■■■), (H) 
where hj = Di{b), k,m eIj and the binomial coefficients are defined as 

. _ n{n - l){n - 2) ■ ■ ■ (n - j + 1) 

This product is associative. 

Definition 5 A formal series 

A = /^D^ + ■ ■ ■ + /o + l-iD-^ + ■ ■ ■ , ken (12) 

is called a formal recursion operator for equation ^ if 

Df{A) = F,oA-AoF,. (13) 

In literature a formal recursion operator is also called a formal symmetry of 
equation ([7j). 

The central result of the Symmetry Approach can be represented by the 
following Theorem, which we attribute to Shabat: 

Theorem 1 // equation ^ has an infinite hierarchy of symmetries of arbi- 
trary high order, then a formal recursion operator exists and its coefficients 
can be found recursively. 

The Theorem states that for integrable equations, i.e. equations possess- 
ing an infinite hierarchy of higher symmetries, one can solve equation ( JT3l) 
and determine recursively the coefficients Im, lm-ii ... of A such that all these 
coefficients will belong to the ring n. The solvability conditions of equation 
ffT^ can be formulated in an elegant form of a canonical sequence of local con- 
servation laws of equation ([71). They provide powerful necessary conditions of 
integrability. These conditions can be used for testing for integrability for a 
given equation or even for a complete description of integrable equations of a 
particular order. A detailed description of the Symmetry Approach including 
the proof of the above theorem and classification results for integrable PDEs 
and systems of PDEs based on the concept of formal recursion operator can 
be found in review papers [IDl [HI [121 fT3] . 
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2.3 Symbolic representation 

Symbolic representation transforms the problems in differential algebra into 
the ones in algebra of symmetric polynomials. This enables us to use powerful 
results from Diophantine equations, algebraic geometry and commutative al- 
gebra. Symbolic representation is widely used in theory of pseudo-differential 
operators. To integrable systems it was first applied by Gel'fand and Dickey 
[21] and further developed in works of Beukers, Sanders and Wang [71 [2S1 12S] • 
Symbolic representation can be viewed as a simplified notation for a Fourier 
transform [29]. 

In order to define the symbolic representation TZ = ©7^2." of the ring (and 
Lie algebra) of differential polynomials TZ = (BTZ^, we first define an isomor- 
phism of the linear spaces TZ^ and TZ^ and then extend it to isomorphisms 
of the differential ring and Lie algebra equipping TZ with the multiplication, 
derivation and Lie bracket. 

Symbolic transform defines a linear isomorphism between the space 7^" of 
differential polynomials of degree n and the space C[^i, . . . , ^n]"^" of algebraic 
symmetric polynomials in n variables, where is a permutation group of 
n variables ^i, . . . , Elements of TZ"^ are denoted by 'u"a(^i, . . . , where 
O'i^ij ■ ■ ■ y^n) € . . . , ^n]"^". The isomorphism of linear spaces TZ^ and 7^" 
is uniquely defined by its action on monomials. 

Definition 6 The symbolic form of a differential monomial is defined as 

u,,u,, ■ ■ ■ G 7^" ^ ■ ■ ■ e 7^" 

where (■)^« denotes the average over the group of permutation of n ele- 
ments ^1, . . . 

For example, 

Uk ^- 
uul 

It is easy to see that in TZ the linear operators Du and Xu, cf. ([3]) are 
represented by 



u 



U , UUi 



(66 + 66 + 66 
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With this isomorphism the hnear spaces corresponding to TZ^ have the 
property that the coefficient functions a(^i, . . . of symbols 

are homogeneous symmetric polynomials of degree p. 

One of the advantages of the symbolic representation is that the action 
of the operator D^, cf. ([2]) is very simple. Indeed, let / G 7^" and / i — > 
M"a(^i, . . . ,in) then 

D.M) ^ M"a(ei,...,U(6 + --- + U, 

and thus ^ M"a(ei,---,en)(ei + --- + a)'- 

Let fen^,f^ uM^u . . . , en) and G 7^'", g ^ ^"6(^1, then 

(i) . The product / ■ g has the following symbolic representation: 

f.g ^ M"a(ei,---,a)°w"Kei,---,u 

= M"+™(a(ei,...,UKa+i,---,W))5e . 

ri + m 

This defines the corresponding multiplication o in TZ. Representation of 
differential monomials (Definition [6]) can be deduced from Uk i — > u^i and 
this multiplication rule. 

(ii) . The Lie bracket [f,g], cf. (jlj) is represented by 

[/, g] ^ u"'^"--' (14) 

{na{^l, . . . , C.n-1, + 1- ^n+m-l)b{$,n, • • • , ^n+m-l) — 

1 H~ ' ' ' H~ ^n+m— l)'^(^m7 • • • 7 ^n+m— l)) c€ 

^n + m— 1 

For example, if / G 7?.\ / t— > 'UCt;(ei) and g G 7?.", (7 ^— ^"0(^1, . . . , ^n), then 
^ (cu(ei + ■ ■ ■ + a) - u;(ei) ^(a)) w"a(ei, . . . ,en). 

In particularly, for f = Ui we have u{C,i) = and [mi, g] = for any g ElZ. 
Thus Ml is a symmetry for any evolutionary equation Ut = g. 

Symbolic representation of differential operators (such as the Frechet 
derivative ([5]) and formal series (ITOll ) is motivated by the theory of linear 
pseudo-differential operators in Fourier representation. To operator Dj. ([2]) 
we shall assign a special symbol r] and the following rules of action on sym- 
bols: 

n 

r/(u"a(ei, . . . ,e„)) = M"a(ei, . . .,in)^ij 
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and the composition rule 



7] o M"a(ei, . . . , a) = uMii, ■■■,UiYl + ^) • 

i=i 

The latter corresponds to the Leibnitz rule o / = Dx{f) + /-Dx- Now it 
can be shown that the composition rule fill I) can be represented as following. 
Let we have two operators fD^. and gD'^^ such that / and g have symbols 
{i*a(^i, . . . , ^i) and u^b{^i, . . . , ^j) respectively. Then 

and 

/D^^D:^^^+^(a(ei,...,e«)(r?+ Umii+i,---,ii+jW)s^ ■ (15) 

m=j+l 

Here the symmetrisation is taken over the group of permutation of all i+j 
arguments ^i, . . . the symbol 1] is not included in this set. In particularly, 
it follows from f|T5l) that o t— > r/'"^^. The composition rule f|T5|) is valid 
for both positive and negative exponents. In the case of positive exponents 
it is a polynomial in rj and the result is a Fourier image of a differential 
operator. In the case of negative exponents one can expand the result on rj 
at ^ oo in order to identify it with ( ITTi) . In the symbolic representation 
instead of formal series (fTOj) it is natural to consider a more general object, 
namely formal series of the form 

B = b{v) + uhi^i, T]) + u%i^i, 6, V) + ^'63(6, 6, 6, r?) + ■ ■ • , 

where coefficients b{ri) ^ 0, 6^(^15 ■ ■ ■ ,^k,v) ^^re rational functions of its argu- 
ments (with certain restrictions which will be discussed in the next section). 

The symbolic representation of the Frechet derivative of the element 
/ I — ^ M"a(^i, . . . ,^„) is 

I — ^ nM''"^a(^i, . . . , 77) . 

For example, let F = + Quui, then F h-* uC,f + 3u'^{^i + ^2) and 

F^^r]^ + 6u(^i + r]) . 

It is interesting to notice that the symbol of the Frechet derivative is al- 
ways symmetric with respect to all permutations of arguments, including the 
argument rj. Moreover, the following obvious but useful Proposition holds 
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Proposition 1 A differential operator is a Frechet derivative of an element 
of TZ if and only if its symbol is invariant with respect to all permutations of 
its argument, including the argument rj. 

The symbolic representation lias been extended and proved to be very 
useful in the case of noncommutative differential rings [28] • In the next 
sections symbolic representation will be extended to the cases of many de- 
pendent variables, suitable for study of system of equations and further to 
the cases of non-local and multidimensional equations. 

2.4 Generalisation to several dependent variables 

The definitions and most of the statements formulated in the previous sec- 
tions [2IIHSII] can be easily extended to several dependents, i.e., to systems 
of equations. In this section we will give a brief account of definitions and 
some results concerning two dependent variables. A generalisation for N de- 
pendent variables is straightforward. For details see [TH [121 1131 EOl [3ll 133] . 
Similarly to w-monomials (Section 12.11) . we define w-monomials 

A differential polynomial / in variables uq, fo, ui, fi, . . . is a finite linear com- 
bination of the form 

(a,/3)eA 

where the sum is taken over a finite set 

A = {(a, /?) I «, /? G Z^o> l«l + > 0}. 

It is a differential ring. We again denote it TZ. Derivation Dx-, (cf. ([2])) is 
now replaced by 

For any f eTZ the Frechet derivative /* is a (row) vector operator 

= = Ef«iV (16) 

\fc>0 fc>0 / 

Systems of two evolutionary equations we will write in vector form 

Uf = F(u,ui, . . . ,u„), (17) 
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where = (m^, VkY and F = (Fi, F2Y are vector-columns where Fi, F2 G 7^ 
(the upper index T stands for the transposition). 

Let us introduce an infinite dimensional linear space over C 

C = {{H^,H2Y\H^,H2en}. 
We equip C with a Lie bracket 

[F,G]=F,(G)-G,(F)e£, 
where the Frechet derivative H* for any H G £ is defined as 

jj ( H\ Hi \ 

Thus C has a structure of an infinite dimensional Lie algebra over C. Subal- 
gebra of symmetries of equation (fT7j) is the centraliser Cc{F) of F in £ (cf. 
section I^TTl) . 

Evolutionary system ( fT7|) defines a derivation Dp -.TZ^TZ 

Df^o) = a*„(Fi) + a*^(F2), 

which is also a derivation of the Lie algebra C This derivation we often 
denote as Df. In this notations derivation coincides with and Du+D^ 
with D-n- 

The ring TZ has several gradings. Here we define a monomial degree 
grading 

7^ = 07^fc, n'' = {aen\D^a = ka}. 

keN 

We say that a = o(7^") if a G 0^^^ 7^^ 

Lie algebra C inherits the gradings of TZ. A monomial degree grading 

C = ^ C\ Cf' = {HeC\ DuH = {k + 1)H}, 

and thus [C^, C^] C /^^'^^ is convenient for the definition of approximate 
symmetries. For H G £ we say that H = o(£") if H G 0fc>„'C'^. Approxi- 
mate symmetries of equation (ITTIl of degree n are defined as elements of the 
approximate centraliser 

C2{F) = {GeC\[F,G] = o{r-')}, 

which is a subalgebra of C. 
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The weighted gradation is useful for the study of homogeneous systems. 
We assign weights w = {wu,Wy) with rational entries to the vector variable 
u and define a linear operator 

W = {wuDu + w^D^ + Xu + X^) 

with spectrum Sw = {(p — + — ^)wv + r | p, g, r G Z>o, p + g > 0}. 
Then the linear subspaces £^ in the decomposition 

£ = £^ = {HG£|l^H = /iH} 

satisfy to 

Cy] C . (18) 

Elements of £^ we call w-homogeneous differential polynomial vectors of 
weight /i. 

For example, if the weight vector of variables {u,v) is w = (1/2, 1), then 
¥=( 

\ U2 + S-Uf 1 + VUl — "iu^ui 

is w-homogeneous element of weight 3/2, indeed F G 

If F G £^ is a homogeneous vector, then is a homogeneous derivation 
of weight /i: 

DfCu C Cy-^-fj, , F G 
From f|T8|) it immediately follows 

Proposition 2 Let G = G,y^ + ■ ■ ■ + G^^ , E be a generator of a 
symmetry of a homogeneous equation, then each w-homogeneous component 
Gj^j, is a generator of a symmetry. 

For evolutionary system ffTTl) a recursion operator (a formal recursion 
operator) A can be defined as a differential or pseudo-differential operator 
(or a formal series) 

A = AfeD^ + Ak-iD^,-' + • • ■ , A, G Mat2x2(7^) , 

which satisfies the following operator equation 

Df(A) = F, oA-AoF, (19) 

(compare with Definition [5]). If action of A is well defined on a symmetry 
Gi, i.e. G2 = A(Gi) G C, then G2 is a new symmetry of the evolutionary 



1 
1 



Wu 

w„ 
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system ffT7|) . Starting from a "seed" symmetry Gi, one can build up an 
infinite hierarchy of symmetries G„+i = A"'(Gi), provided that each action 
of A produces an element of L. 

Symbolic representation of the ring 'R. generated by two independent vari- 
ables M, V is quite straightforward. It is a C-linear isomorphism which is suf- 
ficient to define for the monomials. Suppose we have a monomial u°'v^ . Let 
the symbolic representation for the monomial be 'u'"'a(^i, . . . where 
a(.^i, . . . is a symmetrical polynomial (see Definition [6]). Acting by the 
same rule, but reserving a set of variables Ci, C25 • • • (instead of ^1, ,^2, • • •) for 
the symbolic representation of f-monomials we get 1 — > '0'^'6(Ci, . . . , C|/3|)- 
Then 

^ ^H^l/5la(ei, . . . , eH)6(Ci, . . . , C|/3|) • 

Note that the symbol obtained is invariant with respect to the direct product 
of two permutation groups iS^| x S^^y 

To the product of two elements f,g ^TZ with symbols 

/ t-> ?i"{}"'a(^i, . . . , Ci, . . . , Cm) and g u^v^i^i, . . . , ^p, Ci, • • • , Cg) 

corresponds: 

fg ^ n"+^{;-+^((a(6,..., a, Ci, 

b{^n+l, • • • ; ^n+pj Cm+1; • • • ; Cm+g)) cC ) cC , (20) 

where the symmetrisation operation is taken with respect to permutations of 
all arguments ^ and then ( (the symmetrisation can be made in any order). 

If / G 7^ has a symbol / i — > u"'v"^a{^i, . . . , Ci, . . . , (m), then the 
symbolic representation for the derivative D^if) is: 

DM) ^ U^'v^iil + ■ ■ ■ + a + Cl + ■ ■ ■ + Cm)a(6, . . . , ^n, Cl, • • • , Cm). 

To the operator we shall assign a special symbol rj satisfying the 
following composition rule (the Leibnitz rule) 

77 0u"{;™a(^l,...,^„,Cl,...,Cm) 
= M"{)'"a(ei, . . . , a, Cl, • • • , Cm)(ei + ■ ■ ■ + a + Cl + ■ ■ • + Cm + ^?) • 

For f E 71 with symbol / i — > ■u"{)™a(^i, . . . , Ci, ■ . ■ , Cm) the symbohc 
representation of the Frechet derivative ffTUl) is 

f^u I — ' u"~^{)'"na(^i, . . . , ^n-uV^ Ci, • • • , Cm), 
f^y I — ^ M"w"'"^ma(Ci, . . . , Ci, . . . , Cm-i, v) ■ 



16 



3 Integrability of evolutionary equations 

3.1 Study of symmetries of evolutionary equations in 
symbolic representation 

Using the above symbolic representation of the Lie bracket we can study the 
properties of symmetries in great details. 

Theorem 2 Let the right hand side of evolutionary equation ^ has symbolic 
representation 

F^uu;{^i)+ u^aii^i, 6) + ^'02(6, 6, 6) + ■ ■ ■ 

and the degree of polynomial uj{^i) is greater than 1. // 

G ^ uQi^i) + 6) + ^'^2(^1, 6, 6) + ■ ■ ■ (21) 

is a symmetry, then its coefficients can be found recursively 

Ai(ei,6) = §^||y||ai(ei,6) (22) 

^m-l(^l, ^m) = 77777 ^7 ■■■,Cm)0'm-l{Clj ■■■,Cm) 

^ IQ5 •••5 V 

m—2 j m—l 

i=i \ i=j 
-{m - j)am-i- 

'=0 / 5l ^ 

m m 

G^i^i, = u^iY^U - . (24) 



n=l n=l 



Proof: The proof of the Theorem is straightforward (see, for example [29j). 
Using (1141) we can compute the Lie bracket between F and G. When it 
vanishes up to TZ^, we express Ai{^i,C,2) from the result, which leads to 
formula (J22i) . The Lie bracket vanishing up to TZ"^ is equivalent to formula 

□ 



Corollary 1 For the equation stated in Theorem\^ (i). any symmetry has 
a linear part, that is, fi(Ci) 7^ 0; (ii). algebra of symmetries is commutative. 
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Proof: (i). Let us assume that Q{C,i) = 0. Then it follows from ([22]) that 
Ai{C,i,C,2) = 0. Assuming that = for all 1 < < m — 1, we get from 
( !23l) that Am~i = 0. Thus by induction, we get G = 0. 

(ii). Commutator of two symmetries is a symmetry due to the Jacobi 
identity, but it does not contain a linear part. Thus it must vanish. □ 

Theorem [2] states that a symmetry of equation is uniquely determined 
by its linear part (i.e. dispersion). For fixed all coefficients in the 

series (1211) can be found recursively. Theorem [2] does not mean that any 
evolutionary equation has a symmetry. The right hand side of ( |2T]) must 
represent a valid symbol, i.e., an element of TZ. Thus: 

a. all coefficients Am{^i, ■ ■ ■ ,^m+i) must be polynomial, 

b. it should be a finite number of non-vanishing coefficients Am- 

In general, as it follows from (|22l) and (!23l) . the coefficients Ak are rational 
functions - they have denominators G"^. In order to define a symbol of a 
symmetry, these denominators must cancel with appropriate factors in the 
numerators. Thus factorisation properties of polynomials G"^ and are 
crucial for the structure of the symmetry algebra of the equation. 

Proposition 3 (F. Beukers |38] ) For any positive integer m > 2 the poly- 
nomial 

h — l'f_Lf_Lf_Lf'\'" m-l^m m-ltm m-l^m m-lfm 
ftc,m — [C,l + ^2 + ^3 + — ^1 ~ '^2 S2 ~ ^3 ^3 — C4 ^4 , 

where Uf^iCi ^ 0, is irreducible over C. 

Proof: Suppose that hc^m = A- B with A, and B two polynomials of positive 
degree. Then the projective hypersurface S given by hc,m = consists of two 
components S^, given hj A = 0, B = respectively. The intersection 
Tja f]^B consists an infinite number of points, which should be singularities 
of S since 



dhyQ m dA 



dii dii 



— ■ B 



. dB 



= 0. 



Thus it suffices to show that S has finitely many singular points. 

We compute the singular points by setting the partial derivatives of he 
equals to zero, i.e., 

(^1+6 + 6+^4)'"-'- (01^1)"^-'= 

{^i+^2+^3+ur~'-{c3^3r-' = o 
{^i+^2+^3+^4r-'-{cAur-' = o 
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From these equation follows in particular that 

^1 = Cl/Cl, 6 = C2/C2, 6 = C3/C3, ^4 = C4/C4, 

where C^~^ = 1 and Ci/ci + C2/C2 + Cs/ C3 + C4/C4 = 1. For given q, z = 1 ■ ■ ■ 4, 
we get finitely many singular points. □ 

Corollary 2 Polynomials C^l^i, . . . , c/. ( [^^ are irreducible for n > 4. 

Proof: Let u;(0 = a'^^"' H h ao- If < m < 1 then G"^ is a constant 

and therefore irreducible. If m > 2 polynomial has the form 

G-(ei,...,a) = «,nG(™)(ei,...,a) + ^?^ 

where deg{g'^) < m and 

G^"^) (G, . . . , a) = (^1 + ■ ■ ■ + a)" - C C, (25) 

which is irreducible according to Proposition [3l □ 
Theorem 3 The algebra of symmetries of the evolutionary equation 

n 

Ut = ^akUk + fiun-i,...,u) = F , n>2, a„ 7^ (26) 

k=0 

where /(m„_i, ■ . . ,u) 7^ and 

/K_i,...,«)G0 07^- (27) 

m>3 p<r!, 

trivial, i.e., Cti{F) = SpanQ{ui,F}. 
Proof: In symbolic representation 

F I — > uui^i) + u™+'a^(ei, . . . , U+i) + M'"+'a„+i(ei, . . . , ^^+2) + ■ ■ ■ , 

where 07(^1) = «n^r + ■ — H ai^i + a^. The condition fl271) implies m > 3 and 
degf 

Qm('^l; • • • ) 'Cm+l)) < 

A symmetry of fl2Ul) is of the form 



G ^ Mfi(6) + ^2^1(^1, 6) + ^=^^2(6, 6, 6) + u'Asi^i, 6, 6, ^4) + ■ ■ • 

if it exists. We know that its linear part fi(^i) 7^ from Corollary [H It 
follows from (!23|) that ^^(^i, . . . , ^fc+i) = for /c < m and 

^m(Cl) • • • ) Cm+l) = ^^/^ ' |am(Cl) • • • ) Cm+l)- (28) 
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Suppose ^{^) 7^ + P^^iO for any a,P G C. From Corollary [21 we 
have that polynomials ^^+1) and G'^{^i, ...,^rn+i) are irreducible 

and therefore they are co-prime. Since 

deg(G"^(6, 6m+i)) =n> deg(a„(^i, . . . , ^m+i)), 

the right hand side of ( l28l) is a rational function (not a polynomial). Thus 
there are no symmetries under the assumption. 

When = + /3uj{0 for some e C, it follows from ([23]) that 
G = aui + I3F G Spanj.{Mi, F}. □ 

According to Definition [31 equation ([^ is not integrable. In ([27]) condi- 
tion p < is essential. Indeed, equation 

Ut=U2 + U^u\ 

(p = n = 2) is integrable for any m. 

If an evolutionary equation ([?[) with linear part of order 2 or higher has a 
nontrivial symmetry, then any approximate symmetry of degree 3 is amend- 
able to any degree. Thus, if we have infinitely many approximate symmetries 
of degree 3, then we have infinitely many approximate symmetries of arbi- 
trary high degree. 

Theorem 4 Let uj{^i) he a polynomials of degree greater than 1. Assume 
that evolutionary equation ^ with linear terms uuj{^i) has a nontrivial sym- 
metry. Then for an approximate symmetry X]^=i^j) ^ ^/ degree 3, 
there exists a unique H = X]j>i ^ such that H is an approximate 

symmetry of any degree. 

This Theorem is the direct consequence of a more general Theorem 2.3 
in [26] (see also Theorem 2.76 p. 27 [25]) formulated in the context of filtered 
Lie modules. According to Theorem 2.3, in application to an evolutionary 
equation ([7]) with linear terms uuj{C,i), we should require that polynomials 
G^'^(^i) ^m+i) and polynomials G^{^i, ...,^rn+i) defined by have no 
common factors for some m > 1. This is the case for m = 3 as it follows 
from Corollary [21 

The result of Theorem [H confirms the remark made in [36] : 

Another interesting fact regarding the symmetry structure of evo- 
lution equations is that in all known cases the existence of one 
generalised symmetry implies the existence of infinitely many. 
(However, this has not been proved in general.) 
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For systems of equations and for non-evolutionary equations the conjecture 
that the existence of one symmetry implies the existence of (infinitely many) 
others has been disproved. In [7j it has shown that example given in [6] is 
indeed a counterexample to the conjecture (see also [9]). Even a rectified 
conjecture [39] that for N -component equations one needs N symmetries is 
incorrect either. An example of a system of two equations possessing exactly 
two nontrivial symmetries is given in [H]. 

These examples do not contradict to the spirit of our Theorem H] since 
they are based upon the nonexistence of approximate symmetries of degree 
2, which is one of the conditions in the theorem. 

As we have already mentioned above, the existence of a formal recursion 
operator A ffT^ for an evolutionary equation is a necessary condition for the 
existence of an infinite hierarchy of symmetries. A similar, but not equivalent, 
theorem can be stated in the symbolic representation. The difference is in 
the natural ordering. In the standard representation the coefficients Ik are 
ordered due to the power of Z)^ in the formal series A (fT2!) . In the symbolic 
representation the natural ordering is due to the power of symbol u. The 
fact that Ik must be local, i.e. 1^ ^Ti va. the symbolic representation suggests 
the following definition: 

Definition 7 We say that function bm{^i, ■■■,^m,v)j m > 1 is k-local if in 
the expansion as r] ^ oo 

the first k coefficients Pms{^i, ■■■,^m), s = 1, k are symmetric polynomials 
in We say that bm{^i, ■■■,^m, v) ^■^ local if it is k-local for any k. 

Existence of an infinite hierarchy of symmetries implies the existence of 
a formal recursion operator with local coefficients (Proposition 3 in [29]). 
Existence of an infinite hierarchy of approximate symmetries of degree N 
implies that first — 1 coefficients of the formal recursion operator are 
local. The details of the proof of the following Theorem one can find in [29j 
(Proposition 3). 

Theorem 5 Suppose equation Ij^ has an infinite hierarchy of approximate 
symmetries of degree N 

where fii('^i) '^'^^ polynomials of degree mi and mi < m2 < ■ ■ ■ < rui < ■ ■ ■ . 
Then the coefficients (pmi^i, ■■■,^Tn,f])j ^ = 1,...,A^— 1 of the formal recursion 
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operator 

A = r] + v) + ^^02(6, 6, ^) + ■ ■ ■ 

are local. 

In symbolic representation equation (1131) can be solved ^29j in the sense 
that coefficients of a formal series A can be recursively for any evolutionary 
equation (JTj): 

Theorem 6 Let (f){ri) be an arbitrary function and formal series 

A = 0(r/) + uM^uV) + u^M^u 6, V) + u^M^u 6, 6, ^) + ■ ■ ■ 

be a solution of equation ( fl^) . then its coefficients (f)m{^i, ...,^m,v) can be 
found recursively 

2(0(r/ + ei)-0(r/)) . 

M^uV) = ^j^^ aite,^/) 

<Pm{Cl,--,Cm,V) = ^^f. ^ . r( ("^ + 1) (0(^ + ^1 + ■■■+Cm) 

m—1 

+ ^{n(f>n{^l, ■■,^n-l, H ^ ^m, V)(^m~n{Cn, ••, ^m) 

n=l 

m 
l=n+l 

n . 

-{m-n + l)a„_„(^„+i, ..,^m,V + Y ^l)4>n{^l, •■, ^n, V))s^ 1 • 

Z=l ^ 



Existence of formal recursion operator with local coefficients is a necessary 
condition for the existence of an infinite hierarchy of symmetries. It suggests 
the following test for integrability of equations (JTj): 

• Find a ffist few coefficients (f)n{C,i, ...,^n,v) (ffist three nontrivial coeffi- 
cients (j)n were sufficient to analyse in all known cases to us). 

• Expand these coefficients in series of l/?7 

and check that functions ^nsi^i, ■■■,^n) are polynomials (not rational 
functions). 

This test will be extended and used for non-local and non-evolutionary equa- 
tions in sections l4. II and l 4.2[ 
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3.2 Global classification of integrable homogeneous evo- 
lutionary equations 

In this section, we give ultimate global classification of integrable equations 
of the form 

Ut = Un + f{u,--- ,Un-i), n>2 (29) 

where «„ + /(«,■■■ ,m„_i) is a A-homogeneous differential polynomial and 
A > 0. We give a complete description of integrable equations for all n. 

Theorem H] implies that if equation ( |29l) possesses one higher symmetry 
and infinitely many approximate symmetries of degree 3, then it possesses 
infinitely many approximate symmetries of any degree. Therefore to clas- 
sify integrable equations fl^ it suffices to classify equations, which possess 
infinitely many degree 3 approximate symmetries and then impose the condi- 
tion of existence of at least one exact symmetry. The classification has been 
done in the case of A-homogeneous equations with A > 0. In the case A > 
see the details in while in the case A = the details can be found in 

m 

Now we sketch the results for the case A > without the detailed 
proofs. The following statement on factorisation properties of polynomials 
G^'^^^i, . . . , ^n) (l25l) plays an important role in the classification of integrable 
equations: 

Theorem 7 ^(^^(^i, . . . , = where {g^^\g^^^) = 1 for all k < I, 

and t*^^-* is one of the following cases. 



2: 








k 


= 


(mod 2): 




k 


= 3 


(mod 6).' 


66(6+6) 


k 


= 5 


(mod 6).- 


66(6 + 6)(6' + 66 + ei) 


k 


= 1 


(mod 6).' 


66(6 + 6)(6' + 66 + 6')' 


3: 








k 


= 


(mod 2): 


1 


k 


= 1 


(mod 2): 


(6 + 6)(6 + 6)(6 + 6) 



• n>3: 1 
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For n > 3 the statement follows from more general Theorem El For n = 3, 
it has been proven by Beukers and was published in [23 ES] with his kind 
permission. The case n = 2 has a quite remarkable history. In affine co- 
ordinate X = ^i/^2 we have G^''\^i, ^2) = ^2-^^(2;) and the problem is reducing 
to factorisation properties of the Cauchy-Liouville-Mirimanoff polynomials 

Pkix) = (l + xf-x^-l. 

The common factors Pk = x(l + + x + x'^)'^ g^^\x) and their periodicity 
have been established in the joint report of Cauchy and Liouville [ID]. Using 
Diophantine approximation theory Beukers has shown that factors g^'^^g^"^^ 
are co-prime for k ^ m |41j. Beukers also conjectured that factors g^'^^ are 
irreducible over Q. For prime p the irreducibility of g^'^'^ over Q was earlier 
conjectured by Mirimanoff [12]. A remarkable progress towards the proof of 
the Mirimanoff conjecture has been recently achieved in [13]. 
We now consider A-homogeneous equations of the form 

ut = un + f2 + h + --- , fi^n' (30) 

I — > Mao(6) + ^^^2(6,6,6) H , 

where n > 2, A > and the degree of a polynomial aj is n — jX. Note that if 
A is not integer and iX ^ N, then Oj = 0. This reduces the number of relevant 
A to a finite set. 

Let G e 7^ be a nontrivial symmetry of (l30l) . Then it is of the form 



G = Um + g2 + g3-\ , gi&TV 

^ Mo(6) + M'Ai(ei,6) + M'^2(ei,6,e3) + -- - , 

where 2 < m ^ n and the degree of polynomial Aj is m — jX. For all integers 
r > the following formula holds 

r 

Y,['^'^'a„u'-'+'Ar..]=0. (31) 

i=0 

Clearly we have [uao, uAq] = 0. The next equation to be solved is [uao, u^Ai] + 
[u'^aijuAol = 0, which is trivially satisfied if equation flHUj) has no quadratic 
terms: /2 = 0. Let us concentrate on the case /2 7^ 0. In this case, using 
Theorem HI we see that the existence of a symmetry is uniquely determined 
by the existence of its quadratic term [251 126] . 

We now make a very interesting observation. Assume n and q are both 
odd. Let us compute the symmetry of equation (1301) with linear term Ug. Its 
quadratic terms, cf. (131]) . have the following symbolic expression 
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Proposition [7| implies that A < 3 + 2min(s, s'), where s' = (mod 3) and 
s = ^ (mod 3). We see that if expression (l32i) is a polynomial, then it 
defines a symmetry Q = Uq + Q2 + ■ ■ ■ since Q is determined by its quadratic 
term Q2. The evolutionary equations defined by Q has the same symmetries 
as equation (IHUj) . So instead of (IHUj) we may consider the equation given by 
Q. The lowest possible g is 2s + 3 for s = 0, 1, 2. Therefore we only need to 
consider A-homogeneous equations with A < 7) of orders < 7. 

A similar observation can be made for even n > 2. Suppose we have 
found a nontrivial symmetry with quadratic term 

This immediately implies A < 2. Then the quadratic term 2^iy defines a 
symmetry Q starting with U2. Therefore, we only need to find the symmetries 
of equations of order 2 to get the complete classification of symmetries of A- 
homogeneous scalar polynomial equations (with A < 2) starting with an even 
linear term. 

Finally, we analyse the case when equation fl30l) has no quadratic terms. 
If = for z = 1, ■ ■ ■ , j — 1, then we look at the equation [uao, u^^^Aj] + 
[u^~^^aj,uAo] = 0, i.e. 

_ G(")(ei,---,e,+i)a, 
' GW(ei,---,0+i) ■ 

From Proposition [7] and the proof of Theorem [3], we know there are no sym- 
metries for the equation when j > 3, or when j = 2 and n is even. When 
j = 2 and n is odd, it can only have odd order symmetries. In this case one 
can remark that if the equation possesses symmetries for any m then it must 
possess a symmetry of order 3. 

By now, we have proved the following 

Theorem 8 A nontrivial symmetry of a X-homogeneous equation with A > 
is part of a hierarchy starting at order 2, 3, 5 or 7. 

Only an equation with nonzero quadratic or cubic terms can have a nontrivial 
symmetry. For each possible A > 0, we must find a third order symmetry for 
a second order equation, a fifth order symmetry for a third order equation, a 
seventh order symmetry for a fifth order equation with quadratic terms, and 
the thirteenth order symmetry for a seventh order equation with quadratic 
terms. The last case can be easily reduced to the case of fifth order equations 
by determining the quadratic terms of the equation. The details of this final 
computation are described in [H]. 
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Theorem 9 Let A > 0. Suppose that a X-homogeneous polynomial evolu- 
tionary equation 

Ut = Un + f{u, ■ ■ ■ , Un-l), n>2 

possesses nontrivial symmetries. Then it is a symmetry of one of the follow- 
ing equations up to a transformations u i— >• au, o; e C; 
Burgers equation 

Ut = U2 + UUi 

Korteweg— de Vries equation 

= M3 + UUi 

Potential Korteweg— de Vries equation 

Ut=Uz + u\ 

Modified Korteweg— de Vries equation 

Ut = U^ + U^Ui 

Ibragimov-Shabat equation 

Ut = U^ + 3m^U2 + 9uu\ + ?)U^Ui 

Kaup-Kupershmidt equation 

Ut = u^ + IOMU3 + 2huiU2 + 20m^Ui 
Potential Kaup-Kupershmidt equation 

15 2 20 o 

Ut = U5 + IOM1U3 + + yUi 

Sawada-Kotera equation 

Ut = u^ + lOwMa + 10mi'U2 + 20m^mi 
Potential Sawada-Kotera equation 

20 , 

Ut = U5 + IOM1M3 + 

Kupershmidt Equation 

Ut = u^ + huiu^ + hu\ — hu^u-i — 20 UU1U2 — Swf + bu'^ui 

Finally we note that all the consideration in this section can be extended 
to the case when the dependent variable u takes its values in some associa- 
tive non-commutative algebra (such as matrix, operator, Clifford, and group 
algebras). A complete classification for A > homogeneous equations with 
linear leading term in the case of non-commutative free associative algebra 
was carried out in 1281 . 
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4 Classification results for Non-local integrable 
equations 

The perturbative symmetry approach in the symbohc representation allows 
to derive integrability conditions for certain type of non-local equations. In 
this section we consider two types of such equations: the Benjamin-Ono type 
and the Camassa-Holm type [29l [32] . 
Benjamin-Ono equation reads 

ut = H{u2) + 2uui, (33) 

where H denotes the Hilbert transform 

J-oo y-x 

It is well known that the higher symmetries and conservation laws of the 
Benjamin-Ono equation contain nested Hilbert transform and thus an ap- 
propriate extension of the differential ring TZ is required. The construction 
of such extension is similar to the one proposed by Mikhailov and Yamilov 
in [l5] for 2 -|- 1 dimensional equations (see also section El) . 
The second example is Camassa-Holm type equation: 

rrit = cmui + umi, m = u — U2, c G C \ {0}. (34) 

This equation is known to be integrable for c = 2 [16] and for c = 3 [47] . 
Equation ( IMI) is not in the evolutionary form, but if we exclude one of the 
dependent variables (say u) then we obtain a non-local equation 

mj = cmA(mi) + A(m)mi, A = (1 - D^)"^ (35) 

and again the ring extension is required. 

The symbolic representation and the concept of quasi-locality [13] are the 
key ideas in the extension of the symmetry approach to non-local equations. 
In the definitions of all basic objects such as symmetries, formal recursion 
operators, conservation laws etc. we replace the ring of differential polyno- 
mials TZ by an appropriately extended ring. Elements of this extended ring 
we call quasi-local polynomials (see details in the next section). Symbolic 
representation gives us a simple criteria to decide if a given expression be- 
long to the extended ring. In this extended setting Theorem [2] and most of 
the results of Section 13.11 hold if we just replace "local" by "quasi-local" in 
conditions and statements. 



27 



4.1 Benjamin-Ono type equations and Intermediate long 
wave equation 



Let us consider the following sequence of ring extensions: 

where the set H{JZh") is defined as H{1Zh^) = {H{a)]a G TZh"-} and the 
horizontal line denotes the ring closure. Each TZh" is a ring and the index n 
indicates the nesting depth of the operator H: 

1Z}jo C 7?.j/i C Tl-HZ C ■ • ■ C TZh^ G ■ ■ ■ G TZh^ = T^H- 

Elements of TZh", n > 1 we call quasi-local polynomials. The right hand 
side of equation (IH^ . its symmetries and densities of conservation laws are 
quasi-local polynomials. 

We now consider scalar evolutionary equations, whose right hand side is 
a quasi-local polynomial 

Ut = F, Fe Uh. (36) 

For the definition of its symmetry we replace 7?- by 7^ in Definition [2l Actual 
computations in TZh lead to quite cumbersome calculations. On the other 
hand, in the symbolic representation computations simplify drastically and 
results can be neatly formulated. 

In the symbolic representation operator H is represented by isign{rj). 
So the symbolic representation of the ring extensions is obvious. Suppose 
/ G TZho and 

/ ^ u"a(ei,...,Cn)- 

Then 

HU) ^ u"i sign(ei + ■ ■ ■ + a)a(ei, • • • , ^n). 

In the extended ring all the definitions, such as the Frechet derivative. Lie 
bracket and approximate symmetries, are exactly the same as in the local 
case. However, the symbols of elements of the extended ring are symmet- 
ric sign-polynomials instead of symmetric polynomials. For example, the 
symbohc representation of H{un) and H{uH{ui)) is: 

H{un) ^ Mzsign(^i)^", 

H{uH{ui)) ^ -— sign(ei + 6) (^1 sign(6) + 6 sign(6)) • 
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The symbohc representation of the Benjamin-Ono equation flH^ is 

ut = msign(ei)e? + u\^i + 6)- 

Counting the degrees of sign-polynomials we assume that deg( sign(,^i + ■ ■ ■ + 
6)) = 0. 

Theorem [2] for evolutionary equations and symmetries in TZh, in the sym- 
bolic representation remains the same with the only amendment that all sym- 
bols now are sign-polynomials. To introduce a formal recursion operator for 
equation ( l36l) we introduce a notion of asymptotically local functions, which 
generalises the notion of local functions, cf. Definition [71 

Definition 8 A function a„(^i, . . . , r^) is called asymptotically local if the 
coefficients a„p(^i, . . . , ^n) o,nd a„p(^i, . . . , ^n) of its expansion at t] ^ oo: 

oo oo 

P=Sn P=Sn 

are sign-polynomials, i.e. represent elements from the extended ring TZh. 

In the above expansion we take into account siga{ri + — sign(?7) 

as 7] —>■ oo. 

Definition 9 A formal series 

A = ^{ri) + uM^i, v) + u^M^i, 6, v) + u^M^i, 6, ^3,v) + --- 

is called a formal recursion operator of equation ^3M) if it satisfies equation 
l[T3\) and all its coefficients are asymptotically local. 

Without loss of generality function (j){ri) can be chosen as either ^(r/) = rj or 
(pijj) = rj sign(?7). 

As in the local case, we can solve equation ffT^ with respect to coefficients 
of the formal recursion operator and Proposition [U] holds. The generalisation 
of Theorem [5] is straightforward, we just replace "local" by asymptotically 
local. 

For the Benjamin-Ono equation ( l33i) . the first coefficient (pi^i, rj) of cor- 
responding formal recursion operator 

A = T] + T]) + u'^M^U 6, ^) H 

looks like 

A fc \ ■ f \ , 6(sign(6) + sign(77)) , 1 . 

01 Ui, V) = sign r/ + + O — 

Irj rj' 
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and it is asymptotically local. One may easily check asymptotic locality of 
other coefficients (/>2(Ci, 6, 03(^i, 6, 6, •••• 

In this setting we can classify the generalisation of Benjamin-Ono type 
equation. Consider an equation of the form 

Ut = H{U2) + ClUUi + C2H{uUi) + C2,uH{ui) + C4UiH{u) + 

+c^H{uH{ui)) + CeH{u)H{ui), (37) 

where cj are complex constants. The linear term of this equation coincides 
with the linear term of the Benjamin-Ono equation and all possible homo- 
geneous terms are included if we suppose that if is a zero-weighed operator 
W{H{f)) = W{f) and weight of the variable u equals 2. We also take into 
account that = — 1 and the Hilbert-Leibnitz rule 

H{fg) = fH{g) + gH{f) + H{H{f)H{g)). 

The following theorem holds (see the details and proof in 



Theorem 10 Equation of the form (32) possesses an infinite hierarchy of 
higher syn 
au + hH{u 
of the list 



higher symmetries if and only if it is, up to the point transformation u ^ 
au + hH{u), a'^ + b'^ and re-scalings x ^— ax, t ^ jdt, a, b,a, P G C, one 



Ut = H{u2) + D^{^ciu^ + C2uH{u)-^Ci{uf); (38) 

Ut = H{u2) + D,{^ciu^ + ^C2H{u^)-C2uHiu)y, (39) 

Ut = H{u2) + iuui ± H{uui) =F uH{ui) =F '2uiH{u) 

-iH{uH{ui))- (40) 

Ut = H{u2) + H{uui) + UiH{u) ±iH{uH{ui)) 

±iH{u)H{ui). (41) 

The proof of this theorem requires to check the quasi-locality of the ffist 
three coefficients of the corresponding formal recursion operators. Equations 
dHH]), (HOj) and (HTI) can be reduced to the Burgers equation. In the case 
cl + cl 7^ equation (IHUj) can be transformed into the Benjamin-Ono equation 
(j33ll . When c\ + (^ = 0, it is equivalent to the Burgers equation. The explicit 
form of transformations are in [32]. The properties of the Benjamin-Ono 
equation have been studied in [2]. 

Finally we draw our attention to the intermediate long wave equation 

Ut = -6^^ui + 2uui + T(u2), 
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where 5 is a real constant parameter and 



T{u{x)) ■= ^ J coth (^—{x - y)j u{y)dy . 

This equation was derived by Joseph [IH] as the equation describing propa- 
gation of non-hnear waves in a fluid of finite depth. The intermediate long 
wave equation is an intermediate between Benjamin-Ono and the Korteweg- 
de Vries equations in the sense that the limit 5 — > oo yields Benjamin-Ono 
equation, while 5 — gives KdV equation. The intermediate long wave equa- 
tion possesses an infinite hierarchy of higher symmetries and is integrable by 
the inverse scattering method [19]. As in the case of the Benjamin-Ono 
equation, all its higher symmetries contain nested T operator. 

We consider the general non-linear equation of the intermediate long wave 
form with some linear operator T 

Ut = T{u2) + 2uui (42) 

and address the question: for which linear operators T this equation pos- 
sesses an infinite hierarchy of higher symmetries/conservation laws? In [IH] 
Ablowitz et al have shown that if equation (ll2i) possesses infinitely many 
conservation laws then the conditions 



T{uTv + vTu) = {Tu){Tv) -uv, (43) 




{uTv + vTu)dx = 



must be satisfied. 

The perturbative symmetry approach allows to derive conditions for the 
operator T necessary for the existence of an infinite hierarchy of higher sym- 
metries. All the steps are similar to the case of Benjamin-Ono type equation: 

• Extend the differential ring by the operator T exactly in the same way 
as we did with H and define TZr and its symbolic representation. 

• Define higher symmetries in the extended ring TZr- 

• Introduce a formal recursion operator and asymptotic locality of its 
coefficients. 

Without going into the details (see [22] and [SU]) we present the following 
statement: 
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Theorem 11 Assume that the operator T has the symbolic representation 
tfik): 

/oo 
f{k)u{k)e'^''dk 
'OO 

and that f{k) 1, faster than any power of k~^ , as k +00. Then if equa- 
tion ( [7^ possesses a formal recursion operator with first three coefficients 
being asymptotically local then f{k) satisfies the functional equation: 

fix + y) (fix) + f{y)) = f{x)f{y) + 1 (44) 

Formula fl44l) is equivalent to fl43p in the symbolic representation. Its 
general odd solution, smooth on the real line except the origin, is given by 

f{k) = coth{6k) 

which corresponds to the intermediate long wave equation and the limiting 
case 6 — > +00 corresponds the Benjamin-Ono equation 

f{k) = sign(fc). 

These are the only such equations possessing infinitely many conservation 
laws. The only even solution of fl44l) is 

f{k) = const 

leading to the Burgers equation (up to a re-scaling), which has no non-trivial 
conservation laws. 



4.2 Camassa-Holm type equations 

We now apply the perturbative symmetry approach to determine integrable 
cases of the Camassa-Holm type equation flH31) . It contains operator A = 
(1 — D'^)~^ and therefore we extend the differential ring TZ by operator A 
and define the A-extended ring 71 a as we did for TZh in section 14.11 

The symbolic representation of the operator A is A 1— jz:^- Therefore if 

fen with symbol M"a(^i, . . . , ^„) then A(/) has the symbol w" ^_''i^lX-'+ly' " 
All the definitions remain the same as in the local case with the amend- 
ment 71 — > 7Za- 

We consider a formal recursion operator for the equation fl35|) . First of 
all we introduce a linear term to the equation flH^ by the change of variable 
m I— > m + 1 (note that A(l) = 1) and consider the equation: 

rrit = cA(mi) + cmA(mi) + A{m)mi + mi := F. (45) 
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Its symbolic representation reads: 



where 07(^1) = + ^i, ^(^1, 6) = + gf^y The following state- 



ment holds |32j: 

Theorem 12 The first two coefficients of the formal recursion operator 



The case c = 2 corresponds to the Camassa-Holm equation, while the case 
= 3 corresponds to the Degasperis-Processi equation. In fact, one can show 
[SI] that equation fH31) (or equation flH31) ) with c = 2 or c = 3 possesses an 
infinite dimensional algebra of local higher symmetries in variable m even if 
the equation is non-local. 

5 Integrable Boussinesq type equations 

In this section we give a brief account of our results (see details in [SUl EI]) 
on integrable systems of the form: 



where p>g>r>0, a,l3 E C. System (l46l) can be reduced to a single 
second order (in time) non-evolutionary equation of order p: 

Wtt = aWp + (3Wg^t + F{Wr, Wr+l, Wp_i, Wt, Wi^t,---, Wg-i^t) 

in the variable w, such that v = dtw, u = Wj. and Wk denotes d^w. If function 
F does not depend on v,Vi, . . . , Vr-i then one can eliminate v from the second 
equation and rewrite the system in the form 

Utt = aUp + PUt,q + K{U, Ml, M2, Mt, M<,2, Mt.q-l) , (47) 



A = ?7 + m(j)i{^i,r]) + m^02(^i, 6, ??) H 



for the equation [J5\ ) are quasi-local if and only if c = 2 or c = 3. 




Vt = aUp-r + PVq + F{U, Ml, ...,Up^r~l, V, Ml, Mg_i), 



(46) 



where K = Dl{F). 

The famous integrable Boussinesq equation [521 

Utt Uxxxx ~l~ i^U ")xx 



(48) 
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belongs to this class. Recently all integrable equations of the form 

Utt = Uxxx + Ux, Uxxi Ut, Utx) 

have been classified and comprehensively studied in [5^. In particularly, it 
has been shown that equation 

Wtt = Wxxx + "iWxWtx + WxxWt - 3wlWxx- 

is integrable. 

Sixth order {p = 6) integrable equations of the form (147!) can be obtained 
as reductions of the Sato hierarchies corresponding to KP, BKP and CKP 
equations [Ml ES] as well as derived from equations studied by Drinfeld and 
Sokolov [561 [57] . 

A non-evolutionary equation 

utt = K{u, Ux, Uxx, ■ ■ ■ , d!^u, Ut, Utx, utxx, ■ ■ ■ , d'^ut) , (49) 
can always be replaced by a system of two evolutionary equations 

f Ut = V, 

\ vt = K{u, Ux, Uxx, dxU, V, Vx, Vxx, d^v). 

Non-evolutionary equation (l49l) may have other representations in the 
evolutionary form. U K = Dx{G), then the system of evolutionary equations 

ut = Vx, vt = G . 

also represents fH9l) . 

For example, eliminating variable v from the following systems 
Case 1: 

Ut = V, Vt = Uxxxx + {U^)xx , (50) 

Case 2: 

Ut = Vx, Vt = Uxxx + {U^)x, (51) 

Case 3: 

Ut = Vxx , Vt = Uxx + , (52) 

we obtain the same Boussinesq equation f HHl) on variable u. 

There is a subtle, but important difference among the above represen- 
tations of the Boussinesq equation (HHI) . The system (13UI1 has only a finite 
number of local infinitesimal symmetries (i.e. symmetries whose generators 
can be expressed in terms of u, v and a finite number of their derivatives) 
while equations (15T1) and (|52l) have infinite hierarchies of local symmetries. 
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The reason is simple - infinitesimal symmetries of equation fIFIl) . which de- 
pend on the variable f , cannot be expressed in terms of u and its derivatives, 
since formally v = D~^ut. Equation ( ISTl) possesses only a finite number of 
symmetries that do not depend explicitly on the variable v. 

Homogeneous equations play central role in the theory of integrable equa- 
tions. 

Proposition 4 Let system ( f^^ he w -homogeneous of weight yU, then w = 
Wu + H — r) and 

i. if p is even, p = 2n, n G N, then fi = q = n, 

a. if p is odd, p = 2n + 1, n G N, then P = 0, /i = p/2, q = n + 1. 

System (H6l) we will write in the vector form 

Ui = F, F = {vr,ff. (53) 

Suppose G = {g, h)^ is a generator of a symmetry, then 

[F, G] = ^ Dlih) = Dr{g), D^if) = D^{h). 

It follows from the first equation that D-p{g) belongs to the image of D^. and 
thus the second component of the symmetry generator G can be expressed 
as h = D~^{DF{g)) and the symmetry is completely defined by its first 
component. Substitution into the second equation yields 

DUg) - DlDc^if) = . 

For approximate symmetry of degree n we obviously get 

DUg) - DiDaif) = o{nn . 



Using the above observation for systems ( H6|) we can restrict the action 
of the recursion operator to the first component of symmetries. 

For example, if we represent the Boussinesq equation ( HHi) in the form of 
evolutionary system (1^ . then 



Dl + 2uD^ + 2ui 
and it is easy to verify that a pseudo-differential operator 

A = 



3v + 2viD^^ ADl + 2u + UiD-^ 
A21 3t; + viD-^ 
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where A21 = 4:D^ + lOuDl + 15uiD^ + 9u2 + 4^^ + (2^3 + Auui)D~^, satisfies 
equation ffT^ and therefore is a recursion operator. We can restrict the action 
of the recursion operator on the first component. If gi is the first component 
of a symmetry of the Boussinesq equation (15T|) . then 

g, = ^(S) = {3v + 2v,D-' + {ADl + 2u + u,D-')DtD-') (g,) 

is the first component of the next symmetry in the hierarchy. We call 3? a 
restricted recursion operator. 

A space-shift, generated by ui = {ui,vi)^^, is a symmetry of the Boussi- 
nesq equation fIFIl) . Taking ui as a seed, we can construct an infinite hierarchy 
g3k+i = ^''{ui) of symmetries of weights 3k + 1, k = 0,1,2, . . .. For example, 

(74 = 4f3 + Aviu + Avui = ADx{y2 + vu) . 

We see that (74 is a total derivative and therefore g^ = ^{g^) G 7^ is the 
next symmetry in the hierarchy, etc. The Boussinesq equation itself is not 
a member of this hierarchy. If we take a seed symmetry, corresponding to 
the time-translation {vi,us + 2mmi)^^ we obtain another infinite hierarchy of 
symmetries g3k+2 = ^''{vi),k = 0,1,2,.... The Boussinesq equation does 
not have symmetries of weight 3k, G N. One can show that gsk+i and g3k+2 
are elements of the ring TZ for any /c G N and therefore 3fJ generates two 
infinite hierarchies of symmetries of the Boussinesq equation. Moreover, all 
symmetries from the both hierarchies commute with each other. 

In general, for system ( l53l) a recursion operator is completely determined 
by its two entries An and A12 (the first row of A). The restricted recursion 
operator 3? for system fl53p can be represented as 

3ft = An + Ai2D-'Dt . 

In symbolic representation system fj46l) takes the form 

ut = Kl, 

vt = auCi'" + l3vCl (54) 

It follows from Proposition H] that for w-homogeneous systems (H6ll q = 
p/2 (for even p) and coefficients ajj(^i, . . . , ^j, ^1, . . . , are homogeneous 
polynomials of degree 

degittij) =p+{j - l)r - {i+j - + (j - l)r - ^jp, 

in variables ^1, . . . , ^j, ^1, . . . , Q, they are symmetric in ^1, ... , C,i and in Ci, . . . , Q. 
If deg(ajj) is not a non-negative integer, then = 0. We shall assume that 
the weight Wu > and thus the sum in (15^ is finite due to only a finite 
number of non-zero coefficients an. 
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5.1 Even order equations 

In this section we study homogeneous (with a positive weight Wu > 0) even 
order {p = 2n) equations ( l46i) assuming r = 1. 

In symbolic representation such equations take of the form 

ut = v(i 

vt = wi(^i) +vuj2{Ci) + ^^020(^1, 6) + uvaii{^i,Ci)+ (55) 
+t'^ao2(Ci, C2) + M^a3o(6, 6, 6) + ■ ■ ■ = /, 

where co'i(^i) = o:^i^~^ and ci;2(Ci) = PCi without loss of generality we 
shall represent a in the form 

a = 



4 

and use parameters /i, /3 instead of 

Symmetries of system fl3^ are determined by their first component, which 
in symbolic representation can be written in the form 

= Mfii(ei) + t)^^2(Ci) + ^^20(6, 6) + uvAni^i, Ci) + 

+52^02(0, C2) + m'^30(6, 6, 6) + ■ ■ ■ , 

where f2i(^i), f22(Ci); ^ij are polynomials. For fixed i7i(^i), ^2{Ci) the 
coefficients Aij can be found recursively (there is a generalisation of Theorem 
[2] to the case of many dependent variables [31] ) • 

The Frechet derivative F,, of the system (13H|) with r = 1 in the symbolic 
representation has the form 

7] 

f*,u f*,v 

where 

f*,u = i^iiv) + 2^020(^1, V) + va-niv, Ci) + 3u^a3o(^i, 6, ^) H 

f*,v =^^2(^7) + wan (^1,77) + 2mo2(Ci,^) + M^a2i(6, 6, ^) H 

Formal recursion operator can be defined as a 2 x 2 matrix R whose entries 
are formal series 

^11 = 0OO(^) +M0io(^l,^) + ^^^0Ol(Cl,'7) +^^20(6, 6,^) + 

+uv(l)ni^i,Ci,v) +^^Vo2(Ci,C2,^?) H 

^12 ='ipooiv) + 2^^io(^i>^) +^i'oi{Ci,v) +M^V'2o(6,6,^) + 
+uvi/jii{^i,Cuv) +^^^02(0,(2,??) H , 
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and 



R2i=r] ^ o {Rii^t + o f*,u), R22 = r] ^ o {Ri2,t + R12 o f*,v + Ru ° v)- 
satisfying equation 

Rt = [F*, R] 

and all the coefficients (pij, ipij, i,j = 0, 1, 2, ... of the formal series Ru and 
R12 are local. Approximate recursion operator of degree k can be viewed as 
a truncation of R (terms with ii'^v^ ,i + j > k are omitted or ignored) , so the 
existence of approximate recursion operator is a necessary condition for the 
existence of a formal recursion operator and as well necessary condition for 
the existence of an infinite hierarchy of symmetries of equation ( |55l) . For fixed 
4>oo{v) y i^ooiv) the coefficients of a formal recursion operator can be found 
recursively, similar to the case of one dependent variable (Theorem [U]). The 
property of locality of the coefficients imposes constraints on the coefficients 
of equation (l55l) and eventually leads to the isolation of integrable systems 
(see details in [30]). The latter test had been applied in [30] to the problem 
of classification of even order w-homogeneous integrable systems of the form 
(13^ . The results obtained can be summarised as follows: 



The 4th order equations 

It is easy to see that a homogeneous system ( l46l) with p = 4, r = 1 is linear 
if Wu > S. In the case Wu = S the only possibility is F = ■yu^, 7 7^ which 
leads to a non-integrable equation for any choice of a,f3 and 7. Thus the 
weight Wu can be equal to 2 or 1. 
Case 1: The case of Wu = 2 

The most general nonlinear homogeneous system of equations 0461) with 
p = 4:,r = l,Wu = 2 (correspondingly Wy = 3) is of the form: 

= , (56) 

Vt = aU3 + jJV2 + CiUUi + C2UV, 

where ci, C2 are arbitrary constants and at least one of them is not zero. 

Without loss of generality we need to consider the following three types 
of system fl56|) : 



ut = vi ^ 



Vt = ^ / U3 + PV2 + CiUUi + C2UV, fl ^ {0, ±P} 

= , (58) 

Vt = V2 + CiUUi + C2UV 

(59) 

Vt = -JU3+V2 + C1UU1 + C2UV 
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The first system represents the generic case a ^ {0, — while the other 
two represent the degenerate cases. 

Theorem 13 System (37j possesses two formal recursion operators with 

(pooiv) = V, i^ooiv) = and (pooiv) = 0, i^ooiv) = V 'if and only if p = C2 = 0. 
By re-scalings it can be put in the form 

= , (60) 

Vt = U3 + lUUi 

Systems ( f55l) and ( f5^) are not integrable, they do not possess a formal recur- 
sion operator with 0oo = 0, ipoo = f] unless Ci = C2 = 0. 

System (lUUI) represents the Boussinesq equation (HHI) . which is known to 
be integrable. In the Proposition [T2] and below by re-scalings we mean an 
invertible change of variable of the form: 

uy-^aiu, V ^ a2V, x 1— > as x, t^a^t, G C. 

Case 2: The case of = 1 

The most general homogeneous system of equations fHUl) with p = 4, r = 
1, w„ = 1 is of the form 

Ut = Vl 

Vt = aU3 + [3V2 + C1UU2 + C2U\ + C^UiV + C4UV1 + Csf ^ (61) 
+Cqu'^Ui + C-jU^V + CgW"^ 

where Cj, i = 1 . . . 8 are arbitrary constants and we assume that at least one 
of the coefficients Ci, . . . , C5 is not zero. Without loss of generality we consider 
the following three types of the system (16T|) : 

Ut = Vl 

Vt = ^ ~^ M3 + (iV2 + C1UU2 + C2U\ + C'iUiV + C/iUVi + c^v'^ (62) 

where /i ^ {0, ±/9}, fi, P E C and 
= t;i 

= ^^2 + C1UM2 + C2ul + C3U1V + C4UV1 + Csf ^ (63) 
+Cqu'^Ui + C-jU^V + CgM*^ 

Mf = Vl 

Vt = -|'U3 + V2 + C1UU2 + C2U1 + CsMif + C4Mf 1 (64) 
-l-Csi;^ + CqU^Ui + CrM^W + c%u^ 



39 



Theorem 14 System ^6^) possesses two formal recursion operators with 
000 (^) = V, ^ooiv) = and (poo{r]) = 0, ijjooir]) = r] if and only if (up to 
re-scalings) it is one of the list 



Ut = Vi 
Vt = U3 + u\ 

Ut = Vi 

Vt = Us + 2uiv + 2v?ui 

Ut = Vi 

Ut = Us + 2uiV + Auvi — Gu'^Ui 

Ut = Ui 



(65) 
(66) 
(67) 
(68) 



Ut = Us + Auu2 + Suf — V + Qu ui + u 

Ut = Ui 

Vt = aus + ^2 + 4:auu2 + Sauf + Uiv + 2uvi — f ^ (69) 
+6au'^ui + u^v + au^, a ^ —\ 

System [6&^) possesses a formal recursion operator with 4>oo{v) = 0; i^ooiv) = 
rj if and only if (up to re-scalings) it is one of the list 



Ut = Ui 

Ut = U2 + 2uvi 
Ut = Ui 

Ut = U2 — uf + 2UiV — 
Ut = Ui 

Ut = U2 — 2uu2 — 2u\ + 2uiV + Quvi — V2u^ui 



(70) 
(71) 
(72) 



System possesses two formal recursion operators with (pooiv) = Vt '0oo(^) = 
and tpooiv) = 0; V^oo(^) = V '^f ^'^^ ^^^^2/ ^/ (w re-scalings) it is 

Ut = Ui 

Ut = -jUs + U2 - UU2 - jul + UiV + 2uvi - v'^ (73) 
— Im^-Ui + u'^v — jU"^. 

Equation f l65l) is a potential version of the Boussinesq equation. Equations 
( l66l) and ( 1671) are known to be integrable. Corresponding Lax representations 
and references can be found in [1]. 

Equations (lUSj) and fl69p can be mapped into linear equations 

Wtt = W4 , Wtt = aWi + W2t 

respectively by the Cole-Hopf transformation u = {logw)x- 
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Equation flTUl) can be reduced to the Burgers equation with time inde- 
pendent forcing 

Ut = U2 + 2uui + Wi , Wt = 

by invertible transformation v = w + Ui + u"^ . The latter can be hnearised 
by the Cole-Hopf transformation. 

Equation (ITTl) can be reduced to the system 

Ut = U2 + Wi , Wt = -w'^ 

by a simple invertible the change of the variable v = w + ui. System ( 17T1) 
provides an example of an equation that possesses neither higher symmetries 
nor a recursion operator. However, a formal recursion operator does exist and 
therefore it is in the list of the Proposition. Its integration can be reduced 
to the integration of a linear nonhomogeneous heat equation with a source 
term of a special form. 

By a simple shift of the variable v = w + ui + 2m^, system (1721) can be 
transformed in the form 

Ut = U2 + 4:UUi + wi , Wt = 2Dx{uw) . (74) 

System (1721) possess an infinite hierarchy of symmetries of all orders generated 
by a recursion operator 

3? = -2m + 2uiD-^ + DtD-^ 

starting from the seed Ui. 

Equation flTSj) is a particular case of flUUj) corresponding to the exceptional 
case a = — |. 

The 6th order equations 

Homogeneous 6th order {p = 6) equations (H6|) with non-zero quadratic terms 
correspond to Wu < 5. We restrict ourself with the case Wu > 0. The weights 
3, 4, 5 do not lead to integrable equations: 

Proposition 5 For weights Wu = 3,4,5 there are no equations possessing a 
formal recursion operator. 

Case 1: The case of = 2 

The most general homogeneous system (H6l) with p = 6,r = l,Wu = 2 
can be written as 

Ut = Vi 

Vt = aU5 + I3V3 + Dx{ciUU2 + C2U\ + CsM^) (75) 
+CzUVi + C^VUi, 
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where a,(3,Ci,i = 1,...,4 are arbitrary constants and we assume that at 
least one of Ci, . . . , C4 is not zero. 

Without loss of generality we consider the following three types of the 
above system: 

Ut = Vi 

Vt = ^^-^U5 + (3vs + D^{ciUU2 + C2UI + CgU^) (76) 



Ut = Vi 

Vt = V3 + Dx{CiUU2 + C2U{ + C^U^) + C'^UVi + C4VU1 



(77) 



Ut = Vi 

Vt = -^Ub + V3 + D^{ciUU2 + C2M1 + CsM^) (7J 

+C3UVI + C4VUI 



Theorem 15 If system ( (yq j possesses two formal recursion operators with 
(t>oo{v) = ^) i^ooiv) = '^''T'd 000 (^7) = 0, 'ipooiv) = V then, up to re-scalings, it 
is one of the list 



Ut = Vi 

vt = 2u5 + V3 + D^{2uu2 + ul + ^u^ 

Ut = Vi 

Vt = 1^5 + ^3 + D^{UU2 + UV + \u^) 



Ut = Vi 

Vt = + Vs + D^{2UU2 + + 2uv + |u^) 



(79) 
(80) 
(81) 



If system ([77|j possesses a formal recursion operator with (j)oo{ri) = 0, ipoo{ri) = 
rj then, up to re-scalings, it is one of the list 

(82) 

Vt = V3 + UVi + UiV 
Ut = Vi 

Vt = V3 + 2uus + AU1U2 — Auiv — 8uvi — 2Au^ui 

System ( [7^ does not possess a formal recursion operator for any non-trivial 
choice of Ci . 

Recursion operators and bi-Hamiltonian structure for these equations can be 
found in [30]. Lax representations can be found in [56l [571 EHJ 1301159] . 
Case 2: The case of w„ = 1 
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Homogeneous systems of equations fHUj) with p = 6,r = l,Wu = I can be 
written in the form: 

{Ut = Vi 
vt = au5 + + ciul + C2U1U3 + C3UU4 + C4U2V + c^uiVi ,g^^ 
+CqUV2 + C^v'^ + CgMi + CgUUiU2 + Ciqu'^Us + CnU^Vi + 
+C12UU1V + Ci^u'^ul + CiiU^U2 + Ci^U^V + CigM'^Mi + CirM^, 

where a, /3 G C, all Ci,i = 1, . . . , 17 are arbitrary constants and at least one 
of ci, . . . C7 is not zero. 

We need to consider the following cases of the system fl8^ : 

Ut = Vi 

"^t = ^ ~/ U5 + I3V3 + ClUl + C2U1U3 + CsUU4 + C4U2V 
< +C^UiVi + CqUV2 + Cjv"^ + C^u\ + +C9M-U1-U2 + CiqU^Uj, (85) 
+Cii'U^f 1 + C12UU1V + CiaM^wf + Ci4U^U2 

+ci^u^v + ciqu^ui + cnu^ 
where ^ {0, ±/3}, yU, /5 G C, and 



Theorem 16 If system possesses two formal recursion operators with 
(pooiv) = '4'ooiv) = 0, and 0oo(^?) = 0, ^^00(^7) = V, up to re-scalmgs, it is 




Ut = Vi 

Vt = V3 + Ciul + C2U1U3 + CsUU4 + C4U2V + CsMif 1 

+C(iUV2 + C7f ^ + C^u\ + CgUUiU2 + CioM^Ms + CnM^f i 
+C12UU1V + CisU^uf + Ci4-U^'U2 + Ci^U^V + CiqU'^Ui + CnU^ 



Ut = Vi 

Vt = -^^5 + f3 + ClUl + C2U1U3 + CsUU4 + C4U2V + CsMif 1 
+C6UV2 + Crv'^ + Csul + CgUUiU2 + CioU^Us + CiiU^Vi 
+C12UU1V + Ci^.U^ul + Ci4U^U2 + Ci^U^V + Ci^u'^Ui + CiyM^ 



(86) 



(87) 
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(89) 
(90) 



one of the equations in the following list 

f Ut = Vi 

\ Vt = 2m5 +V3 + ul + 2uiU3 + 
f Ut = Vi 

\ Vt= |m5 + ^3 + UiU^ + UiVi + \ul 
{ Ut = Vi 

\ Vt= |m5 + f3 + 2-U1M3 + |-U2 + 2Mif i + 
Ut = Vi 

Vt = au5 + f 3 + lOaul + Ibauiu^ + Qauu^ + VU2 

+?>uiVi + ?,uv2 — v'^ + Ibaul + \hau^u^+ (91) 
+QQauuiU2 + 'iuuiv + Sw^f 1 + A^av?u\ 

= Vi 

Vt = U5 + QUU4 + ISMiMs + — -u^ + 15-U^'U3+ (92) 

+ 15ul + QOUU1U2 + A5u^ul + 20u^U2 + 15m^ui + 

If system ^8M) possesses a formal recursion operator with (pQo{ri) = 0, ipoo{v) = 
rj, up to re-scalings, it is one of the list 

(93) 

Vt = V3 + UiVi 



Ut = Vi 

Vt = V3 + SuiVi + 3uv2 + Sw^f 1 



(94) 



Ut = f 1 

Vt = V3 - ul + 2U2V - V 



(95) 



If system [8^ possesses two formal recursion operators with (poo{ri) = rj, ipoo{r]) 
0, and 000 (^) = 0, 'ipooiv) = V: W to re-scalings, then it is 



Ut = Vi 

Vt = -\U5 + f3 - \ul - ^MlM3 - \uUi + VU2 ,gg^ 

+3'Uif 1 + 3-uf 2 — f ^ — ^u\ — ^u^U'i — lbuuiU2 + "iuuiv 

+ 3M^f 1 — ^U^u\ — 5u^U2 + U^V — ^u'^Ui — jU^, 

Recursion operators, bi-Hamiltonian structure and Lax representations of 
equations flHHl) . flHI?|) . flUUl) and flUHj) are discussed in details in [5U] . 

Equations fl^ . fP^j) . flMl) and flMj) can be linearised by a Cole-Hopf type 
transformation [30j. 

Equation ( 195|) has similar property as system ( 17T1) . It can also be reduced 
to a triangular system 

ut = U3 + wi, Wt = -w"^ 
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in variables u and w = v — U2- 



10th order equations 

In this section we present three examples of 10th order integrable non- 
evolutionary equations. 

Proposition 6 The following systems possess infinite hierarchies of higher 
symmetries: 



Ut = Vi 

Vt = ^Ug + V5 + [SuUq + 9uiU5 + ^^2^4 + T"3 + 

+2uiVi + Auv2 + 20u^M4 + 8OUU1U3 + 60uul + 88ulu2 

Ut = f 1 

'"t = -^Mg + ^5 + ^UjUi + |m6M2 + 1^5% + 

-S-Usfi - y'^2V2 - WUiVs - ^MsMi - YM1M2M4 

-'iuju, - 'iulu, + 'iulv, + ^«3«? + '^ulul - f 

Ut = Vi 

Vt = + 2M2M5 + 6M3M4 — 6usv — 22-U2'yi — 30'Uit'2 
-20Mf3 + 96uuiv + 96u'^vi + 120D.^{Au^U2 + Gu'^uj) 
-2D^^{8u^Ui + 32mmiM3 + 13m?M2 + 24:uul) - 3840m%i 



(97) 



(98) 



(99) 



Bi-Hamiltonian structures and recursion operators for equations ( 197|) . ( l98l) 
and ( l99i) as well as the Lax representation for equation (1971) can be found in 
[30] . System (p9|l has been considered independently by authors of [591 EH] 
and [51] , where the corresponding Lax representation has been obtained. The 
Lax representations for equation (lUHl) is still not known. 



5.2 Odd order equations 

In this section we formulate the diagonalisation method and globally classify 
homogeneous (with a positive weight Wu > 0) odd order {p = 2n + 1) equa- 
tions ( l46i) in the same spirit as we did for the scalar homogeneous evolution 
equations in Section [31 For the details, we refer the reader to the recent 
paper [HTj. 

A symmetry of an odd order homogeneous equation (1461) always starts 
with linear terms. A homogeneous symmetry in the symbolic representation 
starts either with uC,^ or with v('^~^'^ . Without loss of generality we have 

s 

Ur = G = u^T + J2Yl ^--J-te' -'0, Ci, Cs-,), m > 1 (100) 

s>2 j=0 
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or 

s 

Ur = G = vCr''' + J2J2''''^'~''^^^^-^(^^^-^^^^^^^-^^^-^^^ m>0 (101) 

s>2 j=0 

Functions Aij{C,i, . . . ■ ■ ■ , Cj) iii (HOOP and f llOip are homogeneous poly- 

nomials in their variables, symmetric with respect to arguments Ci, • • • , and 
(i, . . . These functions can be explicitly determined in the terms of sys- 
tem fl54p . the symbolic representation of system (1460 . from the compatibility 
conditions. 

Let us first concentrate on how to compute the Lie bracket between the 
linear part of system ( 15^ denoted by K^, i.e. 

and any pair of differential polynomials. We know its symbolic representation 
takes simple and elegant form if matrix L is diagonal ^23j. Inspired by this, 
we shall diagonalise matrix L, produce the required formula in new variables 
and then transform back to the original variables. 

Notice that matrix L has two eigenvalues ±77""^ 2. Therefore, there exists 
a linear transformation 



T 

such that 



1 1 



T-^LT = diag(77"+^-r7"+5). 
Let us introduce new variables u and v 



V / \ V 



U + V 



T [ ^ 1 — I r).4-i— r n+^—r 



Equally, we have 



u \ _ rp-i ( u \ _ ^ f u + rj" 



The new variables u and -0 have the same weight, i.e., Wu = Wi, = Wu- 
Without causing a confusion we assign the same symbols ^ and ( for the 
symbolic representation of the ring generated by u, v and their derivatives. 
The exponents of symbols can be half-integer, which corresponds to half- 
differentiation in x-space. 
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Proposition 7 In variables u and v, system [54\ ) takes the form 

where g = (^i + ■ ■ ■ + ^« + Ci + " — I" Cs-i)^^^~^ ai^s-i are defined in terms 
ofai^s-iUi,---,^hCi,---,Cs-i), i = 0,...,s as follows: 

s min{i,l} 

ai,s-i{^i, ■■,^1, Cl, ■■, Cs-i) = ^ ^ C'fCs-i 

i=0 p=max{0,l—s+i} 



where are binomial coefficients defined by C/ = -jT(jZjy,- 

In variables u and v, the linear parts of symmetries fllOOl) and (110 II) of 
system (1461) are also diagonal matrices. We can now derive the symmetry 
conditions for the transformed forms as in Theorem[2]for the scalar case. This 
leads to the explicit recursive relations between symmetry (llOOl) or (llOip and 
system (H6l) . 

With explicit formulas at hand, we can prove the following theorem, cru- 
cial for global classification: 

Theorem 17 Assume homogeneous system ( [7^ (p = 2n + 1) with > 
possesses a symmetry. Suppose there is another system of the same weight 
and of the same form 



ut = vCl 

Vt = + Ek>2 E-=o u^v'-%,,.,{^u .., 6, Cl, - , a. 



2n+l-r I v^fc „,i„,k-iu IC c , ^ ^ (102) 



whose quadratic terms equal to those of ^^), that is, bi 2-i{x,y) = ai,2-j(a;, y), i 
0, 1,2. Then if system ^OE) possesses a symmetry of the same order, then 
equation I^102\} and [J^ l (^re equal and sharing the same symmetry. 

This can be viewed as another version of Theorem H] for the systems 
case: the existence of infinitely many approximate symmetries of low degree 
together with one symmetry implies integrability. 

Now we formulate the classification theorem: 
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Theorem 18 If a homogeneous system with a positive weight Wu > 
of odd order (p = 2n + 1) possesses a hierarchy of infinitely many higher 
symmetries, then it is one of the systems in the following list up to re-scaling 
u ^ au, V ^— s> Pv, t 1^ 7t, X ^— s> 6x, where a, P,j,6 E C: 

( Ut = Vi, 

\ Vt = U2 + 3uVi + VUl — Sm^-ui, 
{ Ut = Vi, 

\vt = {D, + uf'^{u)-v\ n = 1,2,3,.... 

These systems can be rewritten in the form of non-evolutionary equations if 
we introduce a new variable u = Wx'- 

Wtt = Wxxx + 3WxWt,x + WxxWt - 3wlWxx- 

The first equation is known to be integrable [53]. The second equation can be 
brought into linear fu = d"^^^ f by the Cole-Hopf transformation w = log(/). 
Its symmetries are given by 

u^^ = DxiDx + u)"'-\ m = 2,3,... 

and 

These symmetries correspond to the symmetries fr^ = fm and fr^ = ft,m of 
equation fu = /2n+i- 

6 Symmetry structure of (2 + 1)— dimensional 
integrable equations 

This section is devoted to the study of (2 + l)-dimensional integrable equa- 
tions. A famous example is the Kadomtsev-Petviashvili (KP) equation 

Ut = Uxxx + QuUx + SD'^Uyy. 

One of the main obstacles to extend the spectacular classification results 
of (1 + l)-dimensional integrable equations to the (2 + l)-dimensional case 
is that the equations themselves, their higher symmetries and conservation 
laws are non-local, i.e. they contain integral operators or D~^. In 1998, 
Mikhailov and Yamilov, [45j, introduced the concept of quasi-local functions 
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based on the observation that the operators D^^ and Dy^ never appear 
alone but always in pairs like D~^Dy and Dy^D,j. for all known integrable 
equations and their hierarchies of symmetries and conservation laws, what 
enabled them to extend the symmetry approach of testing integrability [12] 
to the (2 + l)-dimensional case. 

The results in this section are based on a recent paper We de- 

velop the symbolic representation method to derive the hierarchies of (2 + 1)- 
dimensional integrable equations from the scalar Lax operators and to study 
their properties globally. We prove that these hierarchies are indeed quasi- 
local as conjectured by Mikhailov and Yamilov in 1998, [I5] . 

6.1 Quasi-local polynomials and Symbolic representa- 
tion 

The basic definitions and notations of the ring of (commutative) differen- 
tial polynomials are similar to in section 12. 1[ The derivatives of dependent 
variable u with respect to its independent variables x and y are denoted by 
Uij = dl.d^yU. For small values of i and j, we sometimes write the indices out 
explicitly, that is u^xy and u instead of U21 and Mqo- 

A differential monomial takes the form Ui-^j^Ui^j^ ■ ■ -Ui^j^. We call n the 
degree of the monomial. Let 7^" denote the set of differential polynomials of 
degree n. The ring of differential polynomials is denoted hy IZ = ©„>i7^". 
It is a differential ring with total x-derivation and ?/-derivation 

i,j>0 '-^ i,j>0 

Let us denote 

e = D-^Dy, Q-^ = Dy'Dx. 

The concept of quasi-local (commutative) polynomials TZq was intro- 
duced in |15] . Its definition is similar to the ring extension of TZh as in section 
14. 1[ We consider a sequence of extensions of TZ as follows: 

7^©o = 7^ , Uqu+i = 7e©n |J 0(7^e'O |J Q-^iUen) , 
where the sets 

e(7^en) = {e(/); / e iZe^}, e-l(7^e.0 = {0-^(/); / G 7^e4 
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and the horizontal hne denotes the ring closure. Each TZ^n is a ring and the 
index n indicates the nesting depth of the operators 6^^. Clearly, we have 

7^00 C 7^.01 C 7^.02 C ■ • • C TZ^n C ■ • • C TZqoc = TZq. 

To define the symbolic representation of TZ, we replace Uij by uC,^ri\ where 
^ and 7] are symbols (comparing to the definitions in section [2731) . 

Definition 10 The symbolic representation of a differential monomial is de- 
fined as 

cr&Sn 

The result of action of operators G^^ on a monomial 

in the symbolic representation is given by 

^fc ^ ^ [Vl^ \- Vk\^^ 

M a(^i,...,^fc,?/i,...,77fc) — — . 

V4i H h 

By induction, we can define the symbolic representation of any element in 
7^0, which is a rational function with its denominator being the products of 
the linear factors. The expression of the denominator uniquely determines 
how the operator 0^^ is nested. For example, the symbol of Ux{Qu)Q^^u G 
7^01 is 

V'^^iT J ■ 

3 42^73 4i^3 6^1 

The symbolic representations of pseudo-differential operators in (2 + 1)- 
dimensional case are similar to for the case of one spatial variable in section 
|35j . However, we assign a special symbol X to the operator D^j. in contrast 
to the symbol Y for the operator Dy. 

The extension of the symbolic representation from one dependent variable 
to several dependent variables is straightforward. We need to assign new 
symbols for each of them such as assigning rj^^^ for u and $,^'^\t]^'^^ for v 
and so on. 



50 



6.2 Lax formulation of (2 + 1)— dimensional integrable 
equations 

We give a short description of construction of (2 + l)-dimensional integrable 
equations from a given scalar Lax operator based on the well-known Sato 
approach. For details on the Sato approach for the (1 + l)-dimensional case, 
see the recent books [621 ES] and related references in them. 

Let H be an m-th order pseudo-differential operator in two spatial vari- 
ables of the form 

H = -Dy + a^D™ + a^_iD^-i + ■ ■ ■ + + a^^D'^ + ■ ■ ■ , m > 0, 

where coefficients are functions of x, y. Let the commutator be the 
bracket on the set of pseudo-differential operators. Thus, the set of pseudo- 
differential operators forms a Lie algebra. For an integer A; < m, we split 
into 

H>k = amD^ + ttm-iD^ ^ + ■ ■ ■ + akD^ 
H<k = H — H>k = —Dy + ttk-iD^r^ + • • • 

This operator algebra decomposes as a direct sum of two subalgebras in 
both commutative and noncommutative cases when k G {0, 1}. Similar to 
the (1 -|- l)-dimensional case, such decompositions are naturally related with 
integrability and lead to admissible scalar Lax operators for the case of (2-1-1) 
dimensions: 

a. A; = : n > 2, 

L = + + m("-3)L)^-3 + . . . + -u(o) - Dy; 

b. k = l: n>2, 

L = D2 + u'^^'-^^D^-^ + + ■ ■ ■ + + D-^u^-^^ - Dy; 

c. k = l: L = + D-^u^~^'> - Dy; 

where w^'^ are functions of two spatial variables x, y. We often use u,v,w, - ■ ■ 
in the examples. For the KP equation, the Lax operator is the case a when 
n = 2, namely, L = D^ + u — Dy. 

Let S = D^ + ao + a_iD^^ + ■ ■ ■ . For any operator L listed in cases a, b 
and c, the relation 

[5, L] := SL-LS = 0, (103) 
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uniquely determines the operator S by taking the integration constants to 
be zeros. Furthermore, we have [S"', L] = for any n E N. For each choice 
of 2, we introduce a different time variable tj and define the Lax equation by 

^ = [5V.,L], (104) 

where k is determined by the operator L as listed in cases a, b and c. 

Theorem 19 For the operator S uniquely determined by I^103\) . the flows 
defined by Lax equations (TO^ commute, i.e., dt-dt-L = dt-dt^L. 



6.3 Lax formulation in symbolic representation 

We put the formalism of section 16.21 into the symbolic form. The strategy is 
to do the calculation as much as possible without symmetrisation and only 
perform the symmetrisation at the last stage to get the uniqueness of the 
symbolic representation since the symmetrisation complicates the calculation 
dramatically. 

Let us assign the symbols ^(*),r/(*) for dependent variable u'^^\ The sym- 
bolic representations of the admissible scalar Lax operators are 

a. k = 0: n>2, 

L = X'^-Y + m(«-2)X'^-2 + u{n~3)xn~3 + . . . + ^(0) . 

b. k = l: n>2, 

L = X''-Y + m(«-i)x"-i + ■ ■ ■ + + ; 

c. k = l: L = -r + M(o)+^(-i)_i_. 

Here we only treat the case a. The study of the cases b and c can be found 
in 

It is convenient to consider formal series in the form 

n-2 

S = X + Y,u^%?{^f\v?,X) (105) 

i=0 

n-2 n-2 
«1=0 42=0 

where n > 2 and are functions of their specific arguments, the superindex 
is G {0, 1, 2, ■ ■ ■ , n — 2} and the subindex is defined by the number of ik 
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in the list of [ii, 12, - ■ ■ , ii]- This imphes that ji = 1 and jjt > 1, A; = 1, 2, ■ ■ ■ . 
For example, when ii = 12, the arguments of function 0^2^^ are 
and X. 



It is easy to check that 



. . . ^-^^i) n^* 

where the polynomial A*"; is defined by 



i=l 1=1 



Proposition 8 For any operator L in case a, if the formal series lilU5\) 
satisfies the relation [S, L] = (cf. Iil03\) ). we have for I > I, 

r=l s=r+l 

The function bi, I > 1, is defined by 



biKkn '^ii '^^ - 

Ni{^^\v^\---,^^\V^\X) 



with 



bu^t\vt^---,d^\vt:\x+e'^' 



■=31 ' 131 ' ^3l-l ' 131-1 ' '^31'' 

h /'c(*2) fe) Ail) in) / ^ 1 

-bi-i[U ,Vh '■■■'Q, .Vi, ,X), l>l 



and the initial function Ci{^i\ ri^i\ X) = ^[^\ 

To construct the hierarchy of the Lax equations we need to expand the 
coefficients of operator (11051) at X — 00 and truncate at the required degree. 
When n > 2, the expansion of Ni{C,i, rji, ^2, V2, ■ ■ ■ , ^1, X)~^ aX X ^ 00 is 
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of the form 

k=0 i=0 

Therefore, if we want to prove that the coefficients of operators (11051) . i.e. 
the functions a;, are quasi- local defined similarly as definition [71 we need to 
show that the functions q can be split into the sum of the image of and 
the image of Dy. It is clear that ai\^i\r)^\x) are quasi-local since we have 
ci{^['\r]f\x) = ^['\ When / > 1, we have 

Proposition 9 The functions ci{i^j^\vj^\ ' i '^j?^ ''^j?^ f^'^^ > vanish 
after substitution 

Ah) ^ _Ai2) Ak-i) r ih) ^ fe) 

^jl ^32 ^jl-l 'Ijl 'Ij2 ■ 

In fact, this proposition does not lead to our intended conclusion that bi 
and thus ai are quasi-local since the objects are rational, not polynomial. For 
example, the expression — representing uQu — {Qu)u satisfies the 

above proposition. However, we can not write uQu — {Qu)u = D,j.fi + Dyf2, 
where both /i and /2 are in TZq. 

Notice that the formulae in Proposition [8] are without symmetrisation. 
Combining these expressions, we can obtain the formulae of high degree 
terms of operator S when dependant variables are commuting. Every term 
is such S is quasi-local. This implies that every term in S*" is quasi-local. 
From Theorem [19] follows 

Theorem 20 The hierarchies of commutative (2 + 1 )-integrable equation 
with scalar Lax operators are quasi-local. 

The above setting up is valid for the noncommutative case except The- 
orem [201 However, the extension of the concept of quasi-locality to the 
noncommutative case is rather complicated. Dx and Dy are the only deriva- 
tions for the commutative differential ring. The extension simply enables us 
to apply and Dy^ on the derivations. We know that the commutators 
are also derivations for a noncommutative associative algebra and we need 
to take them into consideration. There are some further discussions on this 
topic in [35] . 



n>2. 
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Summary and discussion 



In this article we have reviewed some recent developments in the symmetry 
approach in the symbolic representation. In particular we have discussed two 
different methods. One method is based on the study of conditions for the 
existence of a formal recursion operator, another one is based on the explicit 
analysis of approximate symmetries of low degrees. 

The formal recursion method is a very efficient tool for testing the inte- 
grability of a given PDE. This method is based only on the fact of existence 
of an infinite hierarchy of higher symmetries and is not sensitive to possi- 
ble lacunas in the hierarchy of symmetries. It is extendable to wide classes 
of equations, including certain types of non-local equations. The method 
is rather simple and convenient for classification of integrable systems of a 
fixed order. We have illustrated its power in applications to classification of 
integrable generalisations of Boussinesq, Benjamin-Ono and Camassa-Holm 
type equations. 

Explicit analysis of approximate symmetries relies on the structure of 
dispersion relations (linear terms) of systems. We have seen that the struc- 
ture of symmetries of a given equation is parametrised by their dispersion 
relations and the analysis of existence of approximate symmetries is based 
on divisibility properties of special polynomials determined by the dispersion 
laws. Such divisibility properties are often obtained via algebraic geometry 
and number theoretic methods. It gives a detailed information on the struc- 
ture of the hierarchy of higher symmetries. So far it is the only method which 
prove to be suitable for a global classification. Here we mean a classification 
of integrable equations in all orders. In the frame of this method it has been 
demonstrated that any integrable homogeneous evolutionary equation 

Ut = Un + f{Un-l,...,u) U > 2, 

where w{u) > is a symmetry (a member of a hierarchy) of one of the 
equation of order 2, 3 or 5 presented in Theorem 

Systems of equations is a considerably more complicated object. Here 
the only "global" result available is a classification of integrable homoge- 
neous Boussinesq type equations of odd order (Section l5.2p . In the theory 
of integrable systems the description of admissible structures of linear terms 
(the dispersion laws) for equations and their symmetries is important and 
yet unsolved problem. It is so called spectrum or dispersion problem. Let us 
consider system of PDEs 

= y4u„ + F(u„_i, . . . ,u), u= [ui,. . . ,un)'^, n>2 (106) 
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where A is a constant N x N matrix. Its higher symmetry, if it exists, is of 
the form 

= B{m)Um + G{Um-l,. . . ,U), TTL > 2 

where B{m) is a constant matrix commuting with A. Let us denote Ai, . . . , Aat 
the eigenvalues of the matrix A and assume that Ai 7^ 0, then the set 
Sa = {A2/A1, . . . , Atv/Ai} is called the spectrum of system fll06p . Similarly 
defined SB{m) = {fJ'2{nT')/ fJ'i{fn), ■ ■ ■ , f^wi^)/ f^iim)} the spectrum of sym- 
metry, where /ii(m), . . . , /iAr(m) are eigenvalues of the matrix B{m). The 
spectrum is invariant with respect to any re-scaling. In many respects it 
reflects properties of symmetries, conservation laws and solutions of the sys- 
tem. For example, the existence of higher conservation laws for a system of 
two even order equations implies that Sa = {—1} (see [20j). 
Integrable system [57] 

J = (5 — 3^/5)u3 — 2uui + (3 — \^)vui + 2uvi + (1 + V5)vvi, 
\ ff = (5 + 3v^)f3 + (1 — \/E)uui + 2vui + (3 + \/E)uvi — 2vvi, 

has spectrum 5*^ = {A2/A1} = { — 1(7 + 3v^)}. It possesses an infinite 
dimensional algebra of symmetries of orders m = 1,3,7,9 mod 10. The 
ratio of parameters /i2(^)//^i(^) is given by 

f,,{m) (l + exp(gf))"^ 
/ii(m) i + exp(^) 

Symmetry approach in symbolic representation can be applied to the 
study of possible spectra and classification of integrable systems. For systems 
of two homogeneous differential polynomial equations of second order the 
problem has been solved in [23]. Recently we have been working on this 
problem for systems of two equations of odd order. The results obtained will 
be published elsewhere. Here we present two rather non-trivial examples of 
integrable systems, which we believe are new. 

The following system 

' ut = (9 - 5\/3)u5 + {2(9 - 5V3)uu2 + (-12 + 7V3)ul} 
+2(3 - V^)u3V + 2(6 - V^)u2Vi + 2(3 - 2V3)uiV2 
-6(1 + V^)uv3 + {2(33 + 19V^)vv2 + (21 + 12\/3)t;f } 
+|(-12 + 7V3)u'^ui + |(3 - 2^/3){vuui 

+1(24 + 13^/3)^2^1 + 1(36 + 20V3)uvvi - f (45 + 26V^)v^vi, 

< 

vt = (9 + 5\/3)w5 + {2(33 - 19V3)uu2 + (21 - 12V3)Mf } 
-6(1 - V^)u3V + 2(3 + 2^/3)u2Vi + 2(6 + VS)uiV2 
+2(3 + V3)uv3 + {2(9 + 5^/3)vv2 - (12 + 7^/3)vf} 
-|(45 - 26^/3)u'^Ul + 1(36 - 20V3)vuui + |(24 - 13V3)u'^vi 
+1(3 + 2V3){v'^ui + uvvi) - 1(12 + 7Vsyvi 
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possesses an infinite dimensional algebra of higher symmetries with 
^i,{m) (l + cxp(f)r 



Hiim) l + exp(^) 
System 



m = 1,5,7,11 mod 12. 



Ut = — 1^5 — Wvvs — 15^1^2 + IOUU3 + 25uiU2 — 6v Vl 

+6v'^ui + 12uvvi — 12-u^-ui, 
Vt = 15f5 + 30f — 30f3-u — 45f2Mi — 35viU2 — lOvu^ 
—Gv'^vi + Gv'^ui + + 12vuui. 

possesses symmetries of orders m = 1, 5 mod 6 with 



/i2(m) (l + exp(f))- 



//i(m) l + exp(^f^) 
There is a reduction v = to the Kaup-Kupershmidt equation. 
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